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Abstract 

Van Holten's covariant algorithm for deriving conserved quantities is presented, with particular 
attention paid to Runge-Lenz-type vectors. The classical dynamics of isospin-carrying particles is 
reviewed. Physical applications including non-Abelian monopole-type systems in diatoms, intro- 
duced by Moody, Shapere and Wilczek, are considered. Applied to curved space, the formalism of 
van Holten allows us to describe the dynamical symmetries of generalized Kaluza-Klein monopoles. 
The framework is extended to supersymmetry and applied to the SUSY of the monopoles. Yet 
another application concerns the three-dimensional non-commutative oscillator. 
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I. INTRODUCTION 



The knowledge of the symmetries is essential in theoretical and condensed matter physics. 
Indeed, symmetries can be exploited to obtain valuable informations on the motion of a 
classical system or after quantization to generate the energy spectrum algebraically. 

The usual classification provides us with discrete and continuous symmetry transforma- 
tions. The discrete symmetries are described by finite groups while continuous symmetries, 
in which we are especially interested, are described by Lie groups. 

A deep basis for the understanding of global conservation laws in modern physics was 
given by Emmy Noether in 1918 [Noether 1918] . She established that conservation laws 
directly follow from the symmetry properties of a physical system. See also [Trautman 19 67]. 
For instance, the invariance by time translation implies the conservation of the energy; the 
invariance by spatial translation yields the conserved momentum and the invariance under 
rotations provides us with the conserved angular momentum. 

In this thesis, we focus our attention on a novel way of deriving conserved quantities 
which has been put forward recently by van Holten |van Holten 2007] . In this formalism, 
invariants are constructed via Killing tensors which are, indeed, the main ingredients of this 
technique. 

Our main endeavor will be to apply van Holten's covariant recipe to various physical 
systems. 

1. Firstly, we clarify the symmetries associated with isospin-Yang-Mills-Higgs field in- 
teractions. To this end, we review, in the context of Kaluza-Klein theories, the 
classical equations describing the motion of an isospin-carrying particle evolving in 
a non-Abelian background. Our presentation follows that of |Kerner 1968j . who first 
introduced these equations, using a "Kaluza-Klein" approach |Kerner 19 68]. 

Next, we discuss the covariant van Holten formalism we use to investigate the sym- 
metries of systems. We note that the symmetry conditions of the van Holten 
formulation are the same as in the Forgacs-Manton-Jackiw (F-M-J) approach 
Forgacs-Manton 1980, iJackiw-Manton 198"U] to symmetric gauge fields. 



2. Most applications of the van Holten algorithm involve various (Abelian but also non- 
Abelian) monopoles and their symmetries. 
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In detail, for a "naked" Dirac monopole, the angular momentum admits a celebrated 
radial term. It has been proved in turn that no globally defined conserved Runge- 
Lenz vector can exist |Feher 1987] . It has, however, been found before by Mcintosh 
and Cisneros, and by Zwanziger (MICZ) [Mcintosh 19701 Zwanziger 1968| that adding 



a suitable inverse- square potential can remove the obstruction such that the combined 
system can accommodate a conserved Runge-Lenz-type vector. 

The archetype of non-Abelian monopoles corresponds to the one introduced in 1968 by 
Wu and Yang in pure Yang-Mills theory [Wu Yang 1968] . O ne can wonder if a particle 



in the Wu-Yang field admits a Kepler-type dynamical symmetry. Generalizing the trick 
of Mcintosh and Cisneros, and of Zwanziger, we find below the most general scalar 
potential such that the combined system admits a conserved Runge-Lenz vector. This 
result had to be expected, since Wu and Yang monopole is in fact an imbedded Dirac 
monopole. 

Although no monopoles were ever seen in high-energy experiments, monopole-like ef- 
fective fields can arise in Condensed Matter Physics. It has been noted by Moody, 
Shapere and Wilczek, for example, that an effective non-Abelian field arises in a di- 
atomic molecule through Berry's phase due to nuclear motion [Wilczek 1986j . For 
some particular value of a certain parameter, it is just a Wu-Yang monopole field. 
For a full range of the parameter, however, the effective field becomes "truly" non- 
Abelian. Electric charge is not more conserved in this case. The system has still 
spherical symmetry, though, and Moody, Shapere and Wilczek do derive a conserved 
angular momentum - but one which has an "unusual" form. But they confess not hav- 
ing a systematic way to obtain it. This goal has been achieved by Jackiw [Jackiw 1986J 
in the F-M-J framework mentioned above. 

Here, after a short outline of Berry's phase, we re-derive the correct expression for 
the conserved angular momentum Ngome 02/2009] , using van Holten's algorithm. In 



addition, we constructed an "unconventional" conserved charge which reduces to the 
square of the electric charge in the Wu-Yang limit. 

3. The next application of van Holten's approach concerns curved spaces of the Kaluza- 
Klein monopole type [Sorkin 19831 ICross 19831 |Gibbons 04/"l986| |Gibbons 12/1986 
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Mimicking what had been done for the MICZ system, we construct, on curved man- 
ifolds, conserved Runge-Lenz-type vectors along the geodesic motion. To this end, 
using the conservation of the "vertical" component of the momentum, we perform a 
dimensional reduction of our curved manifold. This allows us to find the conditions 
under which the dimensionally reduced manifold admits a Killing tensor field associ- 
ated with a Kepler- type dynamical symmetry Ngome 08/2009] . Our strategy is to lift 
3D expressions to the extended Kaluza-Klein manifold. 

Applied to a generalized Taub-NUT metric, we find the most general external potential 
which can be added such that the combined system exhibits a conserved Runge-Lenz- 
type vector. 

In the multi-center metric case Gibbons 12/1986 1, we show that, under cer- 
tain conditions, a conserved scalar of Runge-Lenz-type does exist for two-centers 
Ngome 08/20 09]. For more than two centers no Runge-Lenz-type invariant does exist. 



4. Supersymmetries arise for fermions in a three-dimensional monopole background 
[D'Hoker 1981 |DeJonghe 19951 |Plyushchay 20001 ILeiva 20031 |Avery 2008| . The 
Hamiltonian of the system then involves an additional spin-orbit coupling term, 
parametrized by the gyromagnetic ratio g. 

Below we construct the (super)invariants using a SUSY extension of the van Holten al- 
gorithm. Our clue here is that the symmetry generators can be enlarged to Grassmann- 



algebra- valued Killing tensors |Ngome 03/2010| . Conserved quantities are obtained 
for certain definite values of the gyromagnetic factor : M = 1 SUSY requires g — 2 
Spector|; a Kepler- type dynamical symmetry only arises, however, for the anomalous 



values g = and g = 4. The latter case has the additional bonus to contain an extra 
"spin" symmetry. 

We find that the two contradictory conditions, namely that of having both super and 
dynamical symmetry, can be conciliated by doubling the number of Grassmann vari- 
ables. The anomalous systems with g = and g = 4 will then become superpartners 
inside a unified Af = 2 SUSY system. 

For a planar fermion in any planar magnetic field, i.e. one perpendicular to the plane, 
an M = 2 SUSY arises without Grassmann variable doubling. 



6 



5. We also construct a three-dimensional non-commutative oscillator with no kinetic 
term, but with a non-conventional momentum-dependent potential such that it admits 
a conserved Runge-Lenz-type vector. The latter is derived by adapting van Holten's 



Our system, with monopole-type non-commutativity has the remarkable property to 
confine particle's motion to bounded trajectories, namely to (arcs of) ellipses. The best 
way to figure the motions followed by the particle is to think of them as generalizations 
of the familiar circular hodographs of the Kepler problem, to which they indeed reduce 
when the noncommutativity is turned off. 



method to a "dual" description in momentum space Ngome 
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II. SYMMETRIES AND CONSERVED QUANTITIES IN A NON-ABELIAN 
FIELD THEORY 



The classical equations governing is o spin- carrying particle motion in a non-Abelian 
background are derived using Kerner's Kaluza-Klein framework. The van Holten covari- 
ant method based on Killing tensors and the Forgdcs-Manton-Jackiw approach based on the 
study of symmetric gauge fields are presented. 

A. The "Kaluza-Klein" framework 

In this section, we deal with Kerner's extension of the Kaluza-Klein (KK) approach to a 
non-Abelian gauge theory |Kerner 1 968J. 

Abelian Kaluza-Klein theory 

First of all, let us recall that electromagnetism can be imbedded into general relativity 
(GR) by adding U(l) local gauge invariance to the theory |Kaluza 19191 IKlein 1926] . See 
also [Einstein 19381 IKerner 1981) . Indeed, let us consider the five-dimensional Einstein- 
Hilbert action given by 



Viewing the 5D manifold as a direct product of a 4D space-time with an unobservable 
space-like loop, and assuming that all components of the metric are independent of the extra 
coordinate, y , we get the most general transformations allowed 




(1) 



where G 5 is the coupling constant and 1Z 5 denotes the 5D scalar curvature. 




(2) 



Putting g 44 = V, the 5D metric tensor reads therefore as 




H,v = 0, ••• ,3. 



(3) 



The transformations (pi) imply that transforms as a gauge vector field, 





(4) 
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The "vertical" translation yields hence a U(l) gauge transformation for the vector field 
so that the theory ^ is locally U(l) gauge invariant. The Kaluza-Klein vector can thus 
be identify with the electromagnetic field. 

Let us now embed the metric ^ into the Einstein-Hilbert action defined in ([!]). We have 

det(g AB ) = det(7^)V = g 4 V , 

and it is also useful to calculate the Christoffel connections. The 5D Ricci scalar TZ 5 is 
expressed in terms of the 4D scalar curvature IZ^, the field strength F^ v = d^A u — d v A^ and 
the scalar field V, 

n b = n A - ^vf^f^ - -jy n ^v- ( 5 ) 
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Substituting this 1Z 5 into the action in ([TJ and integrating with respect to the cyclic variable 
y , we obtain the 4D effective action, 

a = -T7rV/ dx A ^g^V(TZ 4 - \VF^ V ) + [ dx'^nVV , (6) 

167tG 5 J Ma 4 87rG 5 J Mi 

where G5 is the 5D Newton coupling constant. The second integral term in (|6| can be 

dropped since it is a surface term and does not affect therefore the equations of the motion. 

Thus we end up with the following 4D action, 

V / ^V^l^-T^), (7) 

wich involves GR and the Maxwell theory, coupled to an additional scalar field V. 

Let us now study the dynamics of a classical point-like test particle of unit mass in our 
5D space-time. Consider 5D geodesic motion, 

where r denotes the proper time. Using the effective theory ^ in Q, a routine calculation 
yields the equations of the motion, 

d iyj^ dx^ ydV\ dq 
dr ^ dr dr dr 

(9) 

d 2 x^ „ dx v dx x „ dx x q 2 d^V 
dr 2 uX dr dr H x dr 2 V 2 

The first equation in ^ tells us that the "charge" , 



is conserved along the 5D geodesies. The latter can also be viewed as being associated 
with translation, in the "extra" direction, generated by the Killing vector d y . The second 
equation in ^ is a 4D geodesic equation involving in addition to the Lorentz force an 
interaction with the scalar field V. See |Kerner 2 000] for a point of view with V = 1 . See 
also [Kibble 19611 ITrautman 1970] . 



Non-Abelian generalization. 



The non-Abelian extension of the 5D KK approach was given by Kerner in [Kerner 1 968] . 
First, we generalize our previous 5D manifold into a (4+d)-dimensional manifold noted as 
M. = M.^ ® S d . The base .M 4 denotes the usual space-time with coordinates x M , and S d 
represents an unobservable (i-dimensional extra space with the locally geodesic coordinates 
y a , a, b = 4, • • • , (3 + d). For definiteness, we fix d = 3 so that S 3 , viewed as a Lie group, 
is isomorphic to the non-Abelian group SU(2). Moreover, the compact manifold S 3 admits 
the isometry generators, Ej = —i^j(y)db, whose algebra reproduces the SU{2) Lie algebra, 

[S i ,S fc ]=ie , ifc S,, (11) 

and which imply the relation, 

e k (y) d b ^(y) - e 3 {y) d b C k (y) = 4 (i/) • ( 12 ) 

The anti-symmetric tensor e l - k denotes the structure constants of the SU (2) non-Abelian 
gauge group. In the KK approach the 7D diffeomorphism symmetry is broken into 4D 
infinitesimal coordinates transformations augmented with translations along the extra di- 
mensions, 

x» — > x» + 6x» , y a — ► y a + f{x")tf(y) . (13) 



Here the f l (x l/ ) are functions. The 7D generalized metric, invariant under (13), then reads 



, liiv + K ah B a B h v B b n ba 
~9cd=\ , C,D = 0,---,6, (14) 



where K a b is the SU (2) invariant metric and 



s; = 4e b . as) 
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The SU(2) Lie algebra- valued one-form A b ^ here will be identified with the Yang-Mills field. 



transforms indeed as a non-Abelian gauge field. Under ( 13 ) the part K ab of the metric (|3 



is preserved. Using the formula — + £,j£ J k if l due to (13), the off-diagonal components 
9ta> of 9cd change as 

A; = Alix) - dj a + el c A\r = A%{x) - DJ a , (16) 

where 

D,r = d,r-e a bc A b J c (17) 
is the gauge-covariant derivative. The field strength of the potential A^ , 

F% = d,A a u - d„Al - el c A\A% . (18) 

changes in turn as 

F; v = F; v - e\J h Fl v . (19) 
For Abelian groups, the structure constants vanish so that the field strength is invariant and 



(17) reduces to simple derivative. 

This is exactly how an infinitesimal gauge transformation, Sy a = f l {x u )^f(y), acts on 



a non-Abelian gauge field. The result (16) differs from the transformation law of Abelian 
gauge fields by the presence of the term e^ c A b f°. 

We now discuss the reduction of the dynamics starting from the 7D Einstein-Hilbert 
action, 

S 7 ~- 



16nG 



V / dVi/v^^. (20) 

TG7 J Ma 



A tedious calculation provides us with the reduced scalar curvature so that the action (20) 
can be reduced as 



167tG 7 Jm a \ vol{S A ) J S 3 4 J 



(21) 



The 7D Newton constant reads G7 = G?/vol(S ) while 7^.4 and IZ3 are the scalar curvatures 



associated with the metrics 7^ and K ab , respectively. The action (21) describes an Einstein- 
like dynamics plus its coupling to the Yang-Mills fields. Note that the second term in (21) 
is given by the curvature of the extra-space. 

We focus our attention on the dynamics of a classical point-like test particle of unit mass 
in (4 + 3)-dimensional space-time. To this end, we consider the Lagrangian for geodesic 
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motion in total space, 



£ = 9cd 



dx c dx D 
dr dr 



For our metric (14), the Lagrange function (22) corresponds to 

dx^ dx v ,dy a An dx^^,dv h 4h dx^^ 

and we evaluate the associated Euler-Lagrange equations, 

, a = 0, • • • , 3 



(22) 



(23) 



d 


( DC ^ 


\ dC 


dr 
< 




' dx a 


d 


( dC \ 


dC 




K 9 m) 


dy c 



(24) 



0, c = 4,5,6. 



The first equation in (24) yields the motion projected into real 4D space-time, whereas the 
second equation describes the motion in 3D internal space. In details, we have 

dx» 



d r C — 2k. 



a l> 



dy b | Ab dx v 
dr v dr 



£°h^A 



_ „ _ dx^ dx^ 
d a C = Oaj,,,,— — + 2n ab 



dr dr 

Let us now identify the following quantity 



bc " dr ' 

dy a Aa dx^ 
dr " dr 



d a A b u 



dx v 
dr 



(25) 



I, 



dy b 



+ A, 



. dx l 



(26) 



dr u dr 

as the classical isospin variable which describes the motion in (non-Abelian) internal space. 



Next, calculating the remaining terms in (24), we find 

dC 



d 



dir. 



dr\d(g)l dr 



d 



dC 



dr \ d {4g) 



A o • A' A r» v j2 
y^fj, lav ~ u v jre/x j j _ j t ^ fav 



(27) 



dr dr 



dr 2 



Collecting the results (25) and (27), we obtain the equations of motion of an isosp in-carrying 
particle in a curved space plus a Yang-Mills field, 



A^rpfi A <T*M A T* ^ A r T^ 

KM Ju „ a KMJU KMJU ,,o .1. KMJU _ 

+ K» — — + l^F h uv l h ^ = 



dr 2 

dir. 



dr dr 



dr 



_ T a Ab d^ _ n 

dr ^^dr 



(28) 



12 



The first equation in (28 ) describes the motion in 4D real space. Note here the generalized 
Lorentz force 

(29) 



7 ^ b dT 



due to the Yang-Mills field, where the electric charge is replaced by the isospin X a . The 
derivation of the latter is analogous to that of the electric charge in the 5D KK theory, since 
it is also the contraction of the Killing vector field generating "vertical" translations with 
the direction field of the geodesic. 



The second equation in (28) says that the isospin is parallel transported in the internal 



space. Remark that the equations (28) can also be obtained using the 7D geodesic equation 

d 2 x c 



dr 2 



+ r c DE — d ^, c p.i-: 



dr dr 



6. 



(30) 



The equations (28) are known as the Kerner-Wong equations. Indeed, some time after 



Kerner, Wong Wong 1970 obtained the same equations by "dequantizing" the Dirac equa- 



tion. Later Balachandran et al. [Balachandran 1977] also deduced the equations (28) using 
a variational principle. Alternatively, they can be studied using a symplectic approach, 
[Duval 1978UDuval 1982llFeher 1986*] . 
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B. van Holten's covariant Hamiltonian dynamics 



The standard approach to identify the constants of the motion associated with 
the symmetries of a given mechanical system is achieved through Noether's theo- 
rem |Forgacs-Manton 1980 , IJackiw-Manton 1980j . summarized in the next subsection. 



More recently, however, an alternative approach has been put forward by van Holten 
[van Holten 2007] . To present his covariant Hamiltonian dynamics, let us consider a non- 
relativistic charged isospin-carrying particle in three-dimensions with Hamiltonian 

7T 2 

H = — + V{x, X a ), n = p-eA. (31) 

Here p and tt define the canonical and the gauge-covariant momenta, respectively, and V 
is an additional momentum-independent scalar potential. The gauge potential A = A a X a , 
with the internal index a = 1, 2, 3 referring to the non-Abelian su(2) Lie algebra, describes 
a static non-Abelian gauge field. Note that all dynamical variables here are gauge invariant. 

Identifying the su(2) Lie algebra of the non-Abelian variable with R 3 , we consider the 
covariant phase space (x, tt, l) , where the dynamics, 

/ = {/,#}. 

is defined by the covariant Poisson brackets, 

( f \ n f dg df df dg df dg 

^ = Djf d^ ~ d^ D > 9 + 61 F ^d7 k - eabc dI^ T ■ (32) 

The field strength and the gauge covariant derivative read 



(33) 



F jk = djA k - d k Aj - e e abc X a A b j A c k , 
D 3 = d 3 -ee ah XA)^ 
respectively. The commutator of the covariant derivatives is recorded as 

[D i ,D j ] = -e abc X a F^ c . (34) 
It is straightforward to obtain the non- vanishing fundamental Poisson-brackets, 

{x l} 7Tj} = 8 ih {vr,, tt,} = eX a F^ {X\X b } = -e abc X c . (35) 
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Let us remark that from the Jacobi identities we can derive the electromagnetic field equa- 
tion, 

fa, fa, n k }} + fa, fa, 7Tj}} + fa, fa, tt,}} = & = . (36) 



We can now derive the Kerner-Wong equations of motion [cf. (28) 





-el a F* 


dx^ 






dt 2 




dt 


dX a 




'dV 


dt 


— £abc£ b ( 





0. 
0. 



(37) 



To construct the dynamical quantities Q(x, 7r,X) which are conserved along the motion, we 
use the covariant van Holten recipe [van Holten 20 07J. The clue here is to expand constants 
of the motion in powers series of the covariant momenta, 



Q{x, 7T,X) = C(x,I) + Ci{x,l)-Ki + —C i j(x,T)'Ki'Kj H 

Requiring Q to Poisson-commute with the Hamiltonian, 

{Q,U} = 0, 

leads us with the set of constraints, 



DiC = eFFffi + djDjV + e abc l a ' 



dd dV 



o(0) 
oil) 



BC ■ BV 

DiC 3 + Djd = eZ a {F? k C kj + F a jk C ki ) + C ijk D k V + e abc l a ^ — , o(2) 

DiCj k + DjC k i + D k Cij = eZ a (F° l Cij k + F^Ci k i + F kl Cuj) + Ci jk iD\V 

ra dC ijk dV 



+£abcF a 



o(3) 



(38) 



(39) 



(40) 



The series of constraints (40) is a priori infinite since the expansion (38) is also infinite. 



But for conserved quantities which admit finite expansion in covariant momenta, the series 



of constraints (40) can be truncated at a finite order- ( n +1) provided we search for order-n 

= , such that the higher-order constraint 



constants of motion. Hence, we can set Cj 1 ...j n+1 



of (40) becomes the covariant Killing equation, 

D {n C V2 ... tn+l) =0, neW 
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(41) 



It is worth noting that apart from the zeroth-order constants of the motion, i.e., which does 
not depend on the covariant momentum, all order-n invariants are deduced from the system- 



atic method (40) implying rank-n Killing tensors, each Killing tensor solving the equation 



(41). These Killing tensors also represent the higher-order coefficient of the expansion (38) 
and, thus, can generate conserved quantities. The intermediate-order constraints of (40) 
determine the other coefficient-terms of the invariant whereas the zeroth-order equation can 
be interpreted as a consistency condition between the scalar potential and the conserved 
quantity constructed. 

To determine the conserved quantities of the system, we can consider the point of view 



consisting on solving first the equation (41) in order to deduce Killing tensors. This task 
is extremely difficult due to the fact that the gauge covariant derivatives do not commute. 
But, in the particular case where 

Ci 2 ...i n+X {x,1) = Ci 2 ... in+1 (x) , (42) 



we can easily solve (41) for n = 1 . We find the general form of the Killing vectors, 

Q = CiY v ..Y m i Yl "' Ym = ciijx 1 + hi , = -dji , (43) 
where bj and the anti-symmetric denote constant tensors. It is worth mentioning that 



from the Killing vectors (43), one can define the associated Killing tensors of Yano-type, 



CiY 1 -Y m (%) , which satisfy the Killing equation 

D^y^.yJx) + Dfi iYl ...Y m {x) = . (44) 

Note that the Killing- Yano tensors are completely anti-symmetric differential forms in their 
indices. 



To solve (41) for n = 2, even under the condition (42), remains, however, an awkward 



task. The trick is to construct the rank-2 Killing tensors CV, ( X^j clS db symmetrized product 



Gibbons 12/1986 of Yano-type Killing tensors, 

Cij(x) = CiYx-Ym CjY v -Y m + CiY v -Y m Cj Yl -Y m • (45) 

As an illustration, consider the two Killing- Yano tensors, 

CiY = UYin 1 , (with n a constant unit vector) , 

CjY = tjYk X k , 
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(46) 



extracted from (43). The symmetrized product (45) of the Killing- Yano tensors (46) 
provides us with the rank-2 Killing tensor generating Kepler-type dynamical symmetry 
Ngome 03/2010| , 



Cij(x) = 28ij (n-x) — riiXj — UjXi . (47) 
In the following, we focus our attention at given Killing tensors. 



For n = 1, (41) provides us with Killing vectors. For example, we have, for any unit 



vector n, the generator of rotations around the axis n, 

C = nxx (48) 

leading to the conserved angular momentum. 

The generator of space translations along an axis n, 

C = an, a e R , (49) 

implies a conserved quantity which is identified with the "magnetic translations" . 



For n = 2, then (41) yields rank-2 Killing tensors. Similarly, for any unit vector ft, 



= 25ij n ■ x — (riiXj + rijXi) (50) 

is a Killing tensor of rank 2 associated with the conserved Laplace- Runge-Lenz vector. 
The rank-2 Killing tensor implying the conservation of energy reads 

Cy = Sij ■ (51) 

The constant rank-2 Killing tensor generating the conserved Fradkin tensor, associated 
with the three-dimensional SU(3) oscillator symmetry, is 

Cij = , = const . (52) 



• For n > 3, the equation (41) provides us with higher-rank Killing tensors which, in 
general, generate product of already known constants of motion. Thus, no new conserved 
quantities and therefore no new symmetries arise in general from these higher-order Killing 
tensors. 



In this section, we outlined the van Holten procedure (40) to derive the symmetries of 
particle in flat space. This recipe can conveniently be extended to curved space provided the 
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partial derivative is replaced by the metric covariant derivative, di — > Vj . Moreover, the 
van Holten algorithm is practical and efficient to derive linear and higher-order invariants 
in the momenta since the only requirement is to have a Killing tensor corresponding to a 
symmetry transformation. 
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C. The Forgacs-Manton-Jackiw approach 



The invariants of a system can also be sought using the Forgacs-Manton-Jackiw approach 
based on the study of symmetric gauge fields |Forgacs-Manton 1980[ IJackiw-Manton 1 980, 



IJackiw 1980j . Let us first consider indeed, the evolution of a free matter system, in the 
absence of a gauge field. In this case, the dynamics is characterized by several conserved 
quantities associated with spacetime diffeomorphisms. For instance, invariance under tem- 
poral translation generates the conserved energy whereas the space rotational invariance 
generates conserved angular momentum. 

Let us now assume that the matter field interacts with an external gauge field, A a . 
In general, the symmetry of the system is broken by this gauge-matter field interaction. 
However, when A a and the matter field are both invariant under the same three-dimensional 
space infinitesimal diffeomorphisms 

5x a = u a , a = 1,2, 3, (53) 

then, the constants of the motion, namely , wins an extra term and can therefore be 
decomposed into two contributing parts, 

^ ^matter ' ^ gauge ' V / 



The first term on the right hand side of (54), which was the total conserved quantity in the 
absence of gauge field, corresponds, in the presence of an external gauge field, to the matter 
contribution augmented with that of the gauge field-matter interaction into the constant of 
the motion, 

^matter = W X , (55) 

where n a represents the gauge covariant momentum. 



The second term of (|54|), namely C£ auge , is the gauge field's additional contribution restor- 



ing the full symmetry (53) of the system. The Forgacs-Manton-Jackiw approach developped 



here is thus a systematic method to construct the gauge field additional contribution C^ auge 
into the constant of the motion. 

Let us first define the Lie derivative of a tensor dragged along the flow, C, described by 
the vector field u a . For this, we consider the infinitesimal coordinate transformation, 

i a ^^ = i a + ^w a , <Jt<l, (56) 
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FIG. 1: 



associated with the diffeomorphisms generate by the vector flow u a . See (53). The formula 



(56) thus implies the tensor field transformations, 

t/ (x a ) — > v'P {x' a ) , A/3 (x a ) — >• A' p (x ,a ) . 
The Lie derivative along C takes the standard form 



(57) 



< = 5A a = uUdpAp + A^u" 

and can be generalize to a (p, g)-tensor field by 



(58) 



j rpmi,...,m p 



U C ni|l,Bi /J,,n2,...,n„ n-i + " • + J nx,...,n q -x,fi n q ^ 



_rp^,m 2 ,...,m p o mi _ _ rpmx,...,m p -x,vp> , ,m p 
± nx,...,n q M ••• I m,...,n, M 



(59) 



Following Forgacs, Jackiw and Manton, we write the symmetry condition of the gauge field, 
, along the vector flow as 

LuA^DaQ", (60) 
where Q w is a differentiable Lie algebra- valued scalar function. Geometrically, the condition 



(60) refers to the action of infinitesimal automorphism of the principal bundle. Note that 



the effect of the gauge freedom on Ap 



A, 



A a — A a + d a A . 



(61) 



does not affect the symmetry condition (60), but only shifts the differentiable scalar field so 
that 

L„A a = D a Q" with g w = Q w + w^A. (62) 
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Consequently, the gauge potential Ap also remains invariant under the symmetry transfor- 
mation generated by u b . Thus, as expected, we can conclude that the symmetry condition 



defined in (60) is gauge invariant. 

An equivalent way to describe the symmetry condition of a gauge field and therefore 
to obtain the gauge field contribution to the constant of the motion is to express the Lie 
derivative in term of the field strength , 

L U A P = ufFtf + Dp (wMj . (63) 



Injecting this result into (60), it is straightforward to obtain the following equivalent sym- 
metry condition implying the gauge field contribution, discussed by Jackiw |Jackiw 1980] . 

F^ = D p C- auge with C g " auge = u»A»-Q». (64) 

Here, the gauge field contribution to the constant of the motion is a differentiable scalar 



function which can be determined by an integration of the equation (64). 

The physics status of the term Cg auge is now clear. Indeed, it represents the response 
of the external (symmetric) gauge field to a spacetime diffeomorphism. It restores the 
symmetry of the system and appears as a Lie algebra-valued scalar field contribution to the 
constant of motion. The complete constant of motion reads therefore as 

C u = u v ir v + J u a {x)F a p{x) dx 13 . (65) 

Let us remark that identifying the Lie algebra of the SU(2) gauge group with R 3 , the 
usual gauge covariant derivative, which we use in this section, becomes the gauge covariant 



derivative defined as (33 ) in the previous section. The rule is simply to replace the generators 
of the Lie algebra, r a (a = 1, 2, 3) , by the components of the isospin vector, X a . Under this 



transformation, the symmetry condition (64) becomes precisely (with no scalar potential) 



the first-order condition in (40) that a linear, in the covariant momentum, conserved quantity 
has to satisfy. 

Thus, the Forgacs-Jackiw-Manton approach is, in fact, equivalent to the van Holten 
procedure for linear invariants. To generalize the first-cited method to higher-order contants 
of the motion, we require the symmetric gauge field to admitting higher-order Killing tensors. 
Then, as in the case of linear conserved quantities, the invariants can, in that event, be 



splitted into the two contributing parts (54), 



^matter ' Cgauge 
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In that event, the matter plus matter-gauge fields contributions give rise to the term 

C^ att er = ^'"^ ■ ■ ■ , (66) 

where u) fll "' fln denotes the Killing tensor field generating the symmetry. The external gauge 
field carries, however, the contribution Cg auge satisfying the constraints, 

D^ 1 C%£"^ = ^"^ . (67) 



We still here recognize the series of constraints of the van Holten algorithm (40) for particle 



evolving in an external gauge field in the absence of an additional scalar potential. 
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III. ABELIAN MONOPOLES 

Dirac's quantization of magnetic monopole strength is obtained from the associativity of 
operators multiplication. (Dynamical) symmetries of the generalized Taub-NUT metric and 
its multi-center extension are investigated. 

A. Dirac monopole 

The concept of magnetic monopole is one of the most influential idea in modern theoretical 
physics. The hypothesis of particles carrying magnetic charge, g , was first made by Dirac 
[Dirac 19 31 J. who observed that the phase unobservability in quantum mechanics allows 
singularities as sources of magnetic fields, just as point electric monopoles are sources of 
electric fields. These singularities define the celebrated "Dirac string" whose position is not 
detectable. This implies the so-called Dirac quantization condition, 



Consequently, the existence of a single magnetic monopole in the universe would explain the 
quantization of electric charge, for which there is no alternative explanation till this day. 

In work preceding Dirac by over fifty years, Maxwell established the equations describ- 
ing the electromagnetism. A surprising asymmetry inside these Maxwell's equations made 
Poincare and J. J. Thompson to infer that a magnetic charge has to be introduced in the 
theory. The Maxwell equations with this assumption then read 



where p e>m and j e>m denote the electric/magnetic charge and current density, respectively. 
But this introduction responded to a mathematical convenience and had, at that time, no 
physical reality; although, at the same period, P. Curie raised the possibility of the existence 
of free magnetic poles |P. Curie 1894j . 

However, studying the motion of a charged particle in the field of an hypothetic isolated 
magnetic monopole, Poincare [Poincare 189 6J observed that, as the particle is no longer deal 
with central forces, the angular momentum is no longer conserved and the motion is no 



eg = he— , N e Z* . 



(68) 




(69) 
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longer necessarily planar. However, a certain amount of angular momentum resides in the 
magnetic field, and that a total angular momentum does exist, 

J = L-q-, L = x x 7f . (70) 

r 

Here L denotes the mechanical angular momentum and the term ( — qx/r) represents 
the Poincare momentum with q denoting the magnetic pole strength. The total angular 



momentum ( 70 ) is conserved along the motion. 



Later, Wu and Yang Wu Yang 1975| showed that the Dirac string, which was introduced 



as a mathematical artifact, can be totally avoided using two different choices of vector 
potential compatible with the monopole field strength. These two patches read 

f — ?-r (l-cos0) for 0<6<- + 5, 

A r = A e = 0,A,= \ ~ 2 (71) 

^-(l + cos#) for --5<9<ir, 

v r sin V 2 

for any arbitrary 5 in the range < 5 < tt/2 . Each region contains a singularity if we try 

to extend them over the entire region around the monopole as Dirac did, but is regular in 

its restricted hemisphere. In the overlapping region 

tt/2 - 5 < 9 < tt/2 + 5, 

the two patches are related by a gauge transformation, 

A N = A s -V(2gA(x)). 

The latter transformation changes the particle wave functions as 

#(x) — ► exp (2iegA(x))^(x) , 

so that requiring the exponential to be single valued everywhere leads to the Dirac quanti- 



zation condition |Wu Yang 1975| , [cf. 68 



From now on, we discuss the Dirac magnetic monopole without reference to singular 
patches or vector potential Jackiw 12/2002| . To this end, we define the Hamiltonian of the 
monopole system as 



7T 2 r 



n= 2^' ^3= "Pi' ~ A 3 > P3 = ~ lhd 3 , ( 72 ) 
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and the following fundamental commutation rules are posited 

[x\ x J ] = 0, [x\ TTj] = iti8), [n, 7Tj] = ie-e ijk B k . (73) 



Taking into account (72) and (73), we derive the gauge- invariant Lorentz-Heisenberg equa- 
tions specifying the motion of a massive charged particle in the external monopole field 
B, 

x = - [K,x = — 

h m (74) 

jf=-\'H,n\ = (vr x B - B x jf) . 
h L J 2mc V ; 

A priori no constraints on the monopole field B are required in the previous equations of 



the motion. Indeed, equations (74) make sense both when B is source-free, V • B = , or 
not, V • B 7^ . However, when we look the Jacobi identities for the commutators of the 
momenta if , we find 

h 2 - - 

Ujk [rri, [itj, TTfc]] = 2e— V • B , (75) 

which vanishes, as it should, for a source-free magnetic fields, B = V x A . 

In order to obtain the exact form of B , we study now the Lie algebra associated with the 
0(3) symmetry of the monopole system. We first remark that the usual angular momentum 
operator, L = x x n , does not satisfy the o(3) Lie algebra, since we get an obstruction term 
inside of the commutator, 

h -> 

[Li, Lj] = ihe ijk L k + ie-e ijk x k (x ■ B) . (76) 

Following Jackiw |Jackiw 1 980J. we restore the spherical symmetry of the system by adding 
a gauge field contribution, C(x) , into the angular momentum, L , 

J = L + C, (77) 

so that we obtain the modified angular momentum algebra, 

[x , J j J = jfcX 

[iTi, Jj] = ihe^T^k + ie— (x{Bj — 5ij(x ■ B) j — ihdiCj (78) 

h 

[J u Jj] = iht ijk L k + ie-e ijk x k (x ■ B) + ihe ljk x m (e k pl e n mp d n Q) . 
It is now clear that the contribution C (x) restores the standard angular momentum algebra, 
[x\ Jj] = ihe l jk x k , [7Tj, Jj] = ihe ijk TT k , [J h Jj] = ihe ijk J k , (79) 
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provided that the following constraints are verified 
d,C j = - (xiB j - 51 (x ■ B) 



C k + -x k (x-B) =0. (80) 
C k = -x k (x-B)- x> {e k e™d m C l ) ° 

Consequently, the conserved generalized angular momentum along the motion becomes 

J = x x jf [x ■ B)x . (81) 



Moreover, the integrability condition coming from the equations in the left hand side of (80), 

[^,^^' = 0, (82) 

— * 

imposes that the field B satisfies the structural equation, 

xAB 1 - Xi d k B j + 5{ (B k + x m d k B m ) - 5{ + x m diB m ) = , (83) 
which can conveniently be solved with the magnetic monopole field, 

B = , (84) 



where g represents the magnetic charge centered at the origin. In fact, the expression (84) 
is the only spherically symmetric possibility consistent with the Jacobi identity (except in 
the origin), 

V-B = 4irg5 3 (x). (85) 

Note that the obstruction occurs only at the isolated location of the magnetic source, at 
origin, which has to be excluded. 



Following Jackiw Jackiw 12/2002] , the violation, at the origin, of the Jacobi identity 



(85) can be better understood by examining the unitary operator, 



U(a) = exp ( — —a ■ tt ) , (86) 



which according to ( 73 ) implements finite translations by a on i, 

U- X {a)x U{a) = x + a. (87) 
As the momenta operators do not commute according to rt73J), we obtain 1 

U(a)U(b) = exp f-^(x,a,^j U(a + b) . (88) 



1 we use the Baker-Campbell-Hausdorff formula: e A e B — e A + B +U A < B ]+- 
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Here 




represents the flux of the magnetic source through the triangle defined by the three vertices 
located at x , x + a and x + a + b . 




FIG. 2: Magnetic flux through the triangle. 



Using (88) it is straightforward to derive the following expressions, 



(u(d)U{b)^j U(c) = exp ^--^$(f,a,6)^ exp ^-^$(x,a + b,c)j U(a + b + c) 
U(a) (u(b)U(cf) = exp - a, 6, c)^ exp (-^®(x,a,b + c) 



U(a + b + c) 



Combining the two previous formulas, we obtain 



(u{a)U{ty U(c) = exp \^-^}(x, a, b, c) J U(a) (u{b)U(cf) , (89) 



where the first term on the right-hand side of (89) reads 

Here Q(x, a, b, c) can be interpreted as the total magnetic flux, 

Q(x,a,b,c) = J dS ■ B = J dxV ■ B , (91) 
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FIG. 3: Magnetic flux through the tetrahedron. 



emerging out from the tetrahedron constructed with the three vectors a , b , c with one 

vertex at x . See Figure [3] below. 

Positing the following axiom of Quantum mechanics; 

Axiom III.l. : Quantum mechanics realized by linear operators acting on a Hilbert space 
requires operator multiplication to be associative. 



We therefore obtain that the phase factor on the right hand side of (89) has to be unob- 
servable so that the flux is quantized for arbitrary a, b and c, 

exp \-^-Q(x , a, b,c)j = exp (-i2nN) = 1 , N eZ. (92) 

Consequently we obtain 

r _ _ h c _ _ 

dxV-B = 2n—N with V • B = 4ng5 3 (x) ^ , (93) 

which provides us with the Dirac's quantization relation |Dirac 1931] . 

14- iVeZ *- < 94 > 



Note that the equation (93) saves the associativity of operators acting on Hilbert space and 
thus implies the quantization of the magnetic charge. The only requirement here is that 
the magnetic field must be a point source or a set of point sources in order to conserve the 



integrality of the left hand side of ( 93 ) 
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Let us now investigate the classical dynamics of a particle evolving in a magnetic monopole 
field, augmented with a scalar potential V . We inquire, in particular, about Kepler-type 
dynamical symmetries. Our investigations lead us to the well-known Mcintosh- Cisneros- 
Zwanziger system [Mcintosh 19701 Zwanziger 1 968 . 



We start with searching conserved quantities associated with the system. A relevant 
recipe to search for constants of the motion is the van Holten algorithm, presented in Section 



II B[ From now on we fix h — c — 1 and we expand the constant of motions in terms of 
covariant momenta, 



Q(x, 7?) = C(x) + Ci(x)7Ti + ■^Cij(x)TTiTTj + 



(95) 



and we require Q to Poisson-commute with the Hamiltonian of the system. This therefore 
implies to solve the series of constraints, 



CAV = 0, 

d i C = eF ij C j + C ij d j V 

diCj + djCi = e (FikCkj + Fj^Cki) + CijkdkV 



o(0) 
(1) 
(2) 



diCjk + djCki + dkCij — e(FuCijk + FjiCiM + F k iCuj) + C^diY o(3) 



Searching for conserved quantity linear in the momentum, we recall that introducing the 
Killing vector generating space rotations, 



C = n x x , 



(96) 



we directly get the associated generalized angular momentum [see (81) 

7 - - x 
J = x x 7r — eg— . 



(97) 



Considering quadratic conserved quantities, we first obtain that the rank-2 Killing tensor, 



ij ij ' 



(98) 



generates the conserved energy of the system, 



£=^ + V(r) 



(99) 
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On the other hand, inserting into the algorithm the rank-2 Killing tensor generating the 
Kepler-type dynamical symmetry, 



(100) 



we solve the second-order constraint with, 



C = eg 



n x x 



(101) 



Next, we insert (100) and (101) into the first-order constraint of the algorithm and we 



investigate the integrability condition of this equation by requiring the vanishing of the 
commutator, 



[d i ,d j ]C = 



A V{r) 



eV 
2r 2 



0. 



(102) 



Thus, the bracketed quantity must satisfy a Laplace equation so that a Runge-Lenz-type 
vector does exist only for radial effective potential of the form, 



eV 

^(r) = — ^- + -+7 with 0,j e R. 



(103) 



Consequently, the zeroth-order constraint is identically satisfied and the solution of the 
first-order constraint reads, 



C = 



n ■ x 



(104) 



Collecting the results (100), (101) and (104), we get the Runge-Lenz vector conserved along 
the particle's motion, 

-_ 1 — — T* 

(105) 



K = -(t? X J - J X 7f) + 



Note that the presence of the fine-tuned inverse-square term in (103) is mandatory for 



canceling the effect of the monopole. For a "naked" monopole, V = 0, in particular, no 
conserved Runge-Lenz vector does exist [Feher 1986*] . 

Now we can give a complete description of the classical motion of a charged particle in 



the Dirac monopole field, augmented with the potential (103). A MICZ system in what 



follows. Firstly, from the angular momentum (97) we obtain 

x{t) 



J ■ 



-eg 



(106) 



30 



so that the trajectory followed by the particle lies on a cone with axis J and fix opening 
angle 9 defined by 



9 = arccos 



-eg 
J 



(107) 



Secondly, the conservation of the Runge-Lenz vector (105) allows us to construct the new 
conserved vector, 



eg 



such that 



N ■ x{t) — J — eV = const 



(108) 



The result ( 108 ) implies that the particle motion also lies in the oblique plan perpendicular 



to N . Consequently, combining (j 1 06|) and (|108|) the particle motion is viewed to be confined 
onto conic sections 



Gibbons 04/1986] , |Feher 02/2009| . 




FIG. 4: The motion lies on the conic section obtained by intersecting the cone, due to conserved 
angular momentum J, with the oblique plane determined by the additional conserved quantity N . 



The particular form of the conic section depends only on the angle (mod [it]) given by 

J-N 



cos (3 



JN 



(109) 



and which determines the inclination of the oblique plane in comparison to the angular 
momentum vector. 
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We thus obtain the following properties : 

/9e 



For 



7T 

0,--a 



(3 = | - a 



the trajectories lie on < 



7T 7T 

a, — 

2 2J 



ellipses 
parabolae 
hyperbolae . 



(110) 



It is worth noting that the momentum-trajectories called the hodographs are also confined 
to a plane perpendicular to the conserved vector N since 



N ■ 7? (t) = . 



Ill) 



But in some interesting cases, the momentum trajectories can be completely determined. 
For example, in the Kepler problem the 7r-trajectories are known to be (arcs) of circle. In the 
context of non-commutative oscillator mechanics (see later), we prove that the hodographs 
of the MICZ-system lie on (arcs) of ellipses |Ngome 06/2010) . 



Another illustration of using the symmetry ( 105 ) is to derive the energy spectrum from 



the dynamical symmetry Feher 10/1986| Horvathy 1990 Feher 02/2009] . To this end, we 



return to quantum mechanics and consider the vectors J and K defined in (97) and (105), 



respectively, as operators of the Hilbert space satisfying the quantized commutation rela- 
tions, 



[Jj, Jj 



ijk<J j 



K\ [Ki,Kj] =i{2 1 -2H)t ijk J k . 



(112) 



Let us define, on the fixed-energy eigenspace H^fr = S^/ , the rescaled Runge-Lenz operator, 

(2 7 -2£)^ K for £< 7 
K = i K for 8 = 7 • (113) 

(2£-2 7 )~5/? f or £> 7 

We therefore obtain the commutation relations between J and K , 

itijkJ k for £ < 7 
for 8 — 7 



[Jj, Jj] itijkJ 



7 K 



(114) 



-i€ijkJ k for £ > 7 
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Thus, 



£ < 7 o(4) Lie algebra 

for { £ = 7 , if and J generate the < (3) (g) R 3 = e(3) (115) 

£ > 7 o(3, 1) Lie algebra 

For a fixed value of the energy, £ < 7 , we consider the more convenient commuting operators 

A = ^(J + K) and B = ^(J - K) , (116) 

verifying the following relations 

[A, A,] = %e l]k A\ [B h BA = ie ijk B k , [A u BA = . (117) 

Then, the operators A and B extend the manifest o(3) symmetry into a dynamical 
o(3) © o(3) = o(4) Lie algebra. The common eigenvector \1/ of the commuting operators, 
H, A 2 , B 2 satisfies, 



A 2 ^ = a(a + , B 2 $> = 6(6 + = £^ , 

where a and 6 are half-integers. Considering the so far non-negative number, 



n 



x/2 7 - 2£ 



we use the Casimir operators, 
-t 2 



K 



■J 2 + e 2 g 2 -l 



2 7 -2£ 



and J • 



V27 - 2£ ' 



;ii8) 



(119) 



(120) 



to obtain the equalities, 



i(a + l) +6(6 + 1) = ^(eV - 1 + n 2 ), 
i(a + l) -6(6 + 1) = (eflr)n. 



(121) 



Solving the equations (j 12 1 [) provide us with the relations, 

2a + 1 = ±(n + eg) 
2b + l = ±(n-eg) 
a — b = ±eg 
a + 6 + 1 = n . 
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(122) 



Let us recall that from equation (94), the product (eg) is quantized in integers or half- 
integers [in units h = c = 1]. Consequently the first relation in (122) implies that n is 



integer or half-integer depending on the value of (eg) being integer or half-integer. 



We can now derive from (119) the bound-state energy spectrum, 



2tV 



£ n = l- sTg, n = ±eg + 1 , ±eg + 2 ■ ■ ■ , 



(123) 



with the integer value of the degeneracy 



71 



2 - e V = (n - eg) (n + eg) 



(124) 



since n and eg are simultaneously integers or half-integers. 
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B. Kaluza-Klein-type monopoles 



Kaluza-Klein theory is one of the oldest ideas attempting to unify gravitation and gauge 
theory |Kaluza 1919] IKlein 1926] . The physical assumption, in this framework, is that the 
world contains four space-time dimensions, plus an extra cyclic dimension so small that it 
can not be observed. Thus, the ordinary general relativity in five dimensions is considered 
to possess a local U(l) gauge symmetry arising from the isometry transformation of the 
hidden extra dimension. 

Later, Sorkin |Sorkin 198 3], and Gross and Perry [Gross 1983] . introduced the Kaluza- 
Klein monopole which is obtained by imbedding the Taub-NUT gravitational instanton into 
Kaluza-Klein theory. The global stationary metric obtained, 

ds 2 = -dt 2 + f{r) (dr 2 + r 2 (dd 2 + sin 2 6d<j) 2 )) + f~\r) (dx 4 + A^f , 

(125) 

with 9 e [0, 7r] , G [0, 2it] , A$ = Dirac potential , 

has lead to exact solution of the equations of the four- dimensional Euclidean gravity, ap- 
proaching the vacuum solution at spatial infinity. 

In 1986, Gibbons and Manton studied the hidden symmetry of Kaluza-Klein-type met- 



rics and found, in the context of monopole scattering Gibbons 04/1986, Gibbons 12/1986 



that the geodesic motion in the Taub-NUT metric admits a Kepler-type dynamical symme- 



try |Feher 10/198~6j IGibbons 19871 IGordani 19881 ICordani 1990] . (See [Feher 02/2009] for a 
review) . 

A better understanding of such hidden symmetries of Kaluza-Klein-type monopoles 
was achieved by various generalizations |Visinescu 01/1994] |Iwai 05/1994] |Iwai 06/1994 



Visinescu 07/1994, Ivan Holten 1991 IGomtetl IVaman 19961 IGotaescu 19991 IGotaescu 20041 



IKrivonos 20061 |Ballesteros 03/2008] |Ballesteros 10/2008] IKrivonos 20091 IVisinescu 20091 
IKrivonos 2010] INersessianl IVisinescu 2011] |Marquette 2011] . 

More recently, Gibbons and Warnick considered geodesic motion on hyperbolic space 
Gibbons 09/2006] and found a large class of systems admitting such a dynamical symmetry. 



Our aim, in this section, is to present a systematic analysis of Kaluza-Klein-type metrics 
admitting a conserved Runge-Lenz-type conserved quantity. To this end, we consider the 
stationary family of metrics, 



dS 2 = f(x) 5 {j dx i dx j + h(x) (dx 4 + A k dx k ) 2 . 



(126) 
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In these metrics, f(x) and h(x) are real functions and the 1-form A k is the gauge potential 
of a charged Dirac monopole. 

Inspired by Kaluza's hypothesis, as the fourth dimension here is considered to be cyclic, 
we use the conservation of the "vertical" component of the momentum to reduce the 
four-dimensional problem to one in three dimensions, where we have strong candidates 
for the way these symmetries act Ngome 08/2009] . Then, the lifting problem can be 



conveniently solved using the Van Holten technique [see section II B 



Let us first investigate the four-dimensional geodesic motion of a classical point-like test 
scalar particle with unit mass. The Lagrangian of geodesic motion on the 4-manifold en- 



dowed with the metric (126) is 



c - 1 '<*> Tt ir + l < ir + A ^ > 2 " u{1) • (127) 

where we also added an external scalar potential, namely U(x), for later convenience. The 
canonical momenta conjugate to the coordinates (x- 7 , x A ) read as 

4 k (128) 

d jOs . / doc dx \ 

Pi= d{dx^Jd£) =h{x){ l^ + Ak ^ )=q - 

The "vertical" momentum, p& = q , associated with the periodic variable, x 4 , is conserved 
and can be interpreted as conserved electric charge. Thus, we can introduce the covariant 
momentum, 

dx 1 

II, ./VM,, (j , /', " '/-I- (129) 

The geodesic motion on the 4-manifold projects therefore onto the curved 3-manifold with 
metric g%j{x) = f(x) dij , augmented with a scalar potential. The Hamiltonian reads as 

1 q 2 
H=- g lJ (x)Ui Uj + V{x) with V(x) = + U(x) . (130) 

For a particle without spin, the covariant Poisson brackets are given by |SouriauT9 70j 

{B, D} = d k B^-^- d k D + qF kl %L ^ , (131) 
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where F^i = d k Ai — diA k is the monopole field strength. Then, the nonvanishing funda- 
mental brackets are 

{x i ,U j } = 6 i j , {U i ,U j } = qF ij . (132) 

We can now deduce the Hamilton equations yielding the geodesic motion of the scalar 
particle on the 3 -manifold, 

x* = {x\ H) = g ij (x)U j} (133) 

n, = u} = q Fa & - diV + r* n fc & . (134) 



Note that the Lorentz equation ( 134 ) involves also in addition to the monopole and potential 



terms a curvature term, typical for motion in curved space, which is quadratic in the velocity. 



We now inquire about the symmetries of the system. For our investigation, we recall that 
constants of the motion, noted as Q , which are polynomial in the momenta, can be derived 
following van Holten's algorithm |van Holten 2007] . The clue in this technique is to expand 
Q into a power series of the covariant momentum, 



Q = c + c i u i + j l c** miij + i n^-n, + 



(135) 



and to require Q to Poisson-commute with the Hamiltonian augmented with an effective 
potential, H = - II 2 + V{x) . This yields the series of constraints, 

C m d m V(x) = o(0) 
d n C = q F nm C m + C n m d m V(x) o(l) 
ViQ + V l C l = q {F vm Cr + F lm 0^) + C a k d k V(x) o(2) 



(136) 



DiCij + VjCu + ViCij — q [Fi m (7y m + Fj m C u m + Fi m Cj/™) 

+ C^pdmVix) o(3) 



which have to be solved. Here the zeroth-order constraint can be interpreted as a consistency 
condition for the effective potential. It is worth noting that the expansion can be truncated 
at a finite order provided some higher-order constraint reduces to a Killing equation, 



T^{h c i2 - in) = 



(137) 
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where the covariant derivative is constructed with the Levi-Civita connection so that 

■n,c' <)>('' • v j ik c k . (138) 

Then, Ci v ..i p = for all p ^ n and the constant of motion takes the polynomial form, 

P-i 1 

Q = Y,TTl Cll ' tkU ^--- U ^ ( 139 ) 



k=0 



The previously presented van Holten recipe is based on Killing tensors of the 3 -manifold. 
Indeed, the conserved angular momentum is associated with a rank-1 Killing tensor (i.e. 
Killing vector), which generates spatial rotations. Rank- 2 Killing tensors lead to conserved 
quantities quadratic in covariant momenta II 's, etc. Note that Killing tensors has been 
advocated by Carter in the context of the Kerr metric [Carter 19 68] . 

Let us discuss two particular Killing tensors on the 3 -manifold which carries the metric, 

giJ (x) = f(x)5 i3 . (140) 

Our strategy is to find conditions for lifting the Killing tensors, which generate the 
conserved angular momentum and the Runge-Lenz vector of planetary motion in flat space, 
respectively, to the "Kaluza-Klein" 4-space. 

• First, we search for a rank-1 Killing tensor generating ordinary spatial rotations as 

C l = g l3 {x)e\ l n k x l . (141) 

We require Cj to satisfy the Killing equation Cj) = , so that we obtain the following 
theorem |Ngome 08/2009] : 



Theorem III.l. On the curved 3-manifold carrying the metric g%j{x) = f(x) Sij, the rank-1 
tensor 

Cj = 9ij( x ) ^ ki nk ^ 

is a Killing tensor generating spatial rotations around the fixed unit vector n , provided 

xxVf(x))-n = 0. (142) 
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Note that Theorem IIII.ll can be satisfied for some, but not all n's. In the two-center 



metric case, for example, it only holds for n parallel to the axis of the two centers (see the 
next section). 

An important case to consider is when the metric is radial, 



f(x) = f(r) 



(143) 



including the Taub-NUT metrics. In that event, the gradient is parallel to x and (142) 



holds for all n's. Thus, Theorem III. 1 is always satisfied for radial metrics. 



• Next, inspired by the known flat-space expression, we consider the rank-2 Killing tensor 
associated with the Runge-Lenz-type conserved quantity 



Cij = 2 9ij (x) n k x k - g ik (x) nj x k - g jk {x) x k 



(144) 



In order to deduce conditions on the metrics admitting a Kepler-type dynamical symmetry, 
we impose T>uC ^ to vanish. A tedious calculation provides us with 



V {i C jt) =2n k x m ( 9ij {x) rf m + g u {x) Y k m + 9jl (x) T k m ) - mx m d mgjl (x) 

- n,j x m d m gu(x) - m x m d m g ij (x) . 



Let us now calculate each term on the right hand side of (145). We first obtain 



rii x m d m gji(x) = / _1 (x) rii g,ji(x) x m d m f(x) 
rij x m d m g u (x) = / _1 (x) rij g u (x) x m d m f(x) 

ni x m d m g i:j {x) = /^(x) n t g^x) x m d m f(x) , 
and next the curvature terms yield, 

2 9jl (x) n k x m T k m = }-\x) n t 9jl {x) x m d m f(x) + f~\x) n m x m 9jl (x) dj(x) 
+ n k 9j i(x) g nk {x) g im (x) x m d n f{x) , 



(145) 



(146) 



2 9a (x) n k x m T k m = f \x) nj 9il (x) x m d m f{x) + f \x) n m x m g u {x) djf(x) 
+ r\x) n k9il (x) g nk (x) g jm (x) x m d n f(x) , 
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2 g l3 (x) n k x m T k m = f~\x) m 9ij (x) x m d m f(x) + f~\x) n m x m 9ij (x) dj(x) 
+ f' 1 (x)n k g ij (x)g nk (x)gi m (x)x m d n f(x) . 



Inserting (146) and the previous curvature terms into (145), we get 



V^Cji) = f 1 (x) \g(ijdi)f(x) n m x m - g(ijXi) n m d m f(x) 
Requiring V^Cji) = yields the following theorem Ngome 08/20091 : 
Theorem III. 2. On the curved 3-manifold carrying the metric gij{x) = f(x) Sij, the tensor 

Cij = 2 g^x) n k x k - g ik {x) rij x k - g jk (x) n { x k 
is a symmetrical rank-2 Killing tensor, associated with the Runge-Lenz-type vector, provided 



n x (^x x V f(x) j = . 



(147) 



Note that the radial metrics (143) satisfy again the Theorem III. 2 so that, in addition 



to the rotational symmetry, they also admit a Kepler-type dynamical symmetry. 

We can also remark that taking into account the compatibility condition of the metric 
tensor, 

V k g i:j {x) = 0, (148) 

the gij{x) always verifies the order- 2 Killing equation V^g^ = 0. Hence, the metric 
tensor is itself a symmetrical rank-2 Killing tensor and the associated conserved quantity is 
the Hamiltonian [Gibbons 19871 Ivan Holten 2007] . 



Having determined the generators of the symmetry which were previously the object of 
our considerations, we can construct the associated constants of the geodesic motion using 



the algorithm (136). We investigate the radially symmetric generalized TAUB-NUT metric 



so that (126) becomes 



dS 2 = fir) 5ij dx i dx j + h(r) (dx 4 + A k dx k ) 



Then, the Lagrangian (|127|) takes the form 



dx k 



Y-U(r), 



(149) 



(150) 
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where the scalar potential U(r) is necessary to furnish Killing 2-tensors. Respectively 
associated with the cyclic variables x A and time t , the conserved electric charge and the 
energy read 

q 2 



, , k / doc , doc \ 



£ 



IP 



+ 



+ U{r) . 



2/(r) 2h(r) 

Using the relations ( 151[ ), we can rearrange the dimensionally reduced Hamiltonian as 



(151) 



U = - n 2 + /(r) W{r) with W(r) = U{r) + 



T 



+ 



£ 



£ 



(152) 



2" ' ' v v ' 2 //!/•) /!/•) 

which we can now use to derive conserved quantities via the algorithm ( |136 ). 

1) First, we look for conserved angular momentum which is linear in the covariant mo- 

• • • = so 



mentum since the 3-metric now satisfies Theorem III. 1 Hence, C, 



that (136) reduces to 



' C m d m (f(r)W(r)) =0 o(0) 

d n C = qF nm C m o(l) 

^ ViQ + ViQ = . o(2) 
The second- and the first-order constraints yield 

Q = g im (r) e m nk n n x k and C = -qg n k 



(153) 



x 



(154) 



respectively. The zeroth-order consistency condition in (153) is satisfied for an arbitrary 



radial effective potential, providing us with the conserved angular momentum, 



J = x x II — qg 



x 



(155) 



involving the typical monopole term. 

2) Let us now turn to quadratic conserved quantities. In that event, we have 
which implies the series of constraints, 

C m d m (/(r) W{r)) = o(0) 

d n C = qF nm C m + C™d m (f(r)W(r)) o(l) 

V i C l + V l C i = q(F im C l m + F lm C i m ) o(2) 



(156) 



ViCij + VjCu + ViQj = 
41 



o(3) 



Taking = gij(r) as a rank-2 Killing tensor, we deduce from the second-order 



equation of (156) that Cj = 0. As expected, the first-order and the zeroth-order consistency 



relation are both satisfied by any radial effective potential C — f(r) W(r) . The conserved 
energy associated, therefore, read as 

1 



£ = ^IT 2 + fir) W(r) . 



(157) 



Next, we search for a Runge-Lenz-type vector generating the Kepler-type dynamical 



symmetry of the system. Since Theorem HI. 2 is satisfied by the considered radial 3 -metric, 



we have to solve the constraints (156) using the rank-2 Killing tensor 



Cij = 2g ij (r)n k x k - g ik (r)n,jX k - g jk (r)riiX k 



(158) 



inspired by its form in the Kepler problem. We solve the second-order constraint of (156) 
and we get 

qg 



Ci — — gim(r) e j k v? x . 
Next, inserting (467) and ([468 ) into the first-order constraint of (156), we obtain 



(159) 



diC= , (f(r) W(r))' + fl] Xjntxt _ i m w{r) y + <,V 



n 



j ■ 



It is now easy to analyze the integrability condition of the previous equation by requiring 
the vanishing of the commutator, 



[d i ,d j ]C = 



A f(r)W(r) 



q 2 g 2 

2r 2 



0. 



(160) 



Thus, the bracketed quantity must satisfy the Laplace equation. 

The zeroth-order equation is identically satisfied. Consequently, a Runge-Lenz-type con- 
served vector does exist only when the radial effective potential is 



f(r) W{r) 



^ + -+7 with /3, 7 E R. 
Zr r 



(161) 



Equivalent to the result of Gibbons and Warnick jGibbons 09/2006] , the formulas (152) and 
(161) allow us to announce the theorem Ngome 08/20091 : 



Theorem III. 3. For the generalized TAUB-NUT metric (149), the most general potentials 



U (r) permitting the existence of a Runge-Lenz-type conserved vector are given by 



U{r) 



H h 7 



2r^ 



r 



f(r) 2h(r) 



+ £ 



(162) 
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where q and g denote the particle and the monopole charge. And (5 , 7 are free constants 



and £ is the fixed energy [cf. (151)] 



Inserting now (161) into the first-order constraint of (156) provides us with 



which is solved by 



Pi n — — - — k 

c n o n n — rij~ x x n , 



C = - n k x k . 
r 



(163) 



(164) 



Collecting the results (|467|), (|468|) and (|474|) yield the conserved Runge-Lenz-type vector, 

(165) 



x 



K = Ux J + (3-. 

r 



Due to the simultaneous existence of the conserved angular momentum (155) and the con 



served Runge-Lenz vector (165), we obtain a complete description of the motion for gener- 



alized TAUB-NUT metric. Indeed, the motions of the particle are confined to conic sections 
Feher 10/1986] . Our class of metrics, which satisfy Theorem III. 3 includes the following. 



1. The original TAUB-NUT case |Sorkin 1983[ IGross 1983] with vanishing external 
U(r) = 0, 



f(r) 



1 



h(r) 



1 + 



4m 



(166) 



where m is real Feher 10/1986, IGibbons 1987] . We note that the monopole scat- 



tering case corresponds to m — —1/2, see Gibbons 04/1986 Feher 10/1986 



Gibbons 12/1986 . We then obtain for 



7 = q 2 /2 — £ and charge g = ±4m . 
the conserved Runge-Lenz vector, 

K = Ux J -Am(£ -q 2 )- . 



(167) 



(168) 



2. Lee and Lee |Lee 2000] argued that for monopole scattering with independent com- 



ponents of the Higgs expectation values, the geodesic Lagrangian (127) should be 
replaced by L — > L — U(r), where the external potential reads 



U(r) 



1 a 2 



2 4m 
1 + — 

r 



(169) 
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It is now easy to see that this addition merely shifts the value in the brackets in (160) 
by a constant and corresponds to a shift of a^/2 in the energy. Hence, the Laplace 



equation in (160) is still satisfied. So the previously found Runge-Lenz vector (168) is 
still valid. 



3. The metric associated with winding strings [Gibbons 1988] where 

/(r) = l, h(r) = - ' 



(170) 



For charge g = ±1 , we deduce from Theorem HI. 3 



{13 + q 2 ) - r U{r) - 7 + - - f ) = 



2 

so that for the fixed energy, £ = q 2 /2 — 7 + U(r) , the conserved Runge-Lenz vector 
reads as 



K = x x J — q' 



x 



(171) 



4. The extended TAUB-NUT metric [Iwai 05/1994] |Iwai 06/1994] where 

/(r) = »+!!, h(r) = , l r , + i f , , 
r 1 + ar + cr z 



(172) 



with the constants (a, b, c, d) e R. With the choices f/(r) = and charge g = ±1 
requires 



Theorem III. 3 



r /(r) (j 2 q 2 

-r f{r) 8 + — — — - 7 r = /3 = const . 

n(r) 2 2r 



Inserting here (172) yields 



{-a£+ l -dq 2 -p)+r(-b£+ l -cq 2 - 7) = 0, 



which holds when 



f3 = —a£ + -dq 2 and r ) = —b£ + -cq 2 



(173) 



Then, we get the conserved Runge-Lenz vector 



K — II x J — [ a£ dq } - 



x 



(174) 
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5. Considering the oscillator-type metric discussed by Iwai and Katayama Iwai 05/1994 
Iwai 06/1994] , the functions f(r) and h(r) take the form 

ar 4 + br 2 



f(r) = b + ar and h(r) 



(175) 



1 + cr 2 + dr 4 

A direct calculation leads to the following Runge-Lenz-type vector |Marquette 2 010 1 , 

K= (b + ar 2 ) f x J + /3-. (176) 

Which is conserved only for scalar potential of the form 

'1 + cr 2 + dr 4 * 



U(r) 



ar 4 + 6r 2 



(177) 



Let us just conclude by outlining that the five examples treated above are shown to be 
particular cases deduced from the general expression (165), see Ngome 08/20091 . See also 
Igata 20 10 1 for applications on the Kerr metric. The case of SUSY of the Kerr metric is 
investigated in |Galajinsky] . 
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C. Mult i- center metrics 



The multi-center metrics family in which we are interested in this section are known to 
be Euclidean vacuum solutions of the Einstein equations, with self-dual curvature. The 
multi-center metrics can also be viewed as an extension of the TAUB-NUT metrics studied 
in the previous section |Gibbons~9 87j. 

Let us begin by considering a scalar particle moving in the Gibbons-Hawking space 



Gibbons 01/1979 , which generalizes the TAUB-NUT space. The Lagrangian function as- 



sociated with this dynamical system, as in (127), is given by 



But here the functions /(f) obey the following "self-dual" ANSATZ |Gibbons 01/1979] , 

V/ = ±VxA. (178) 



Hence, /(f) is an harmonic function, 



A/(f) = 0. 



(179) 



The most general solution of (179) is given by 

with (f , mi) e R N+l 



N 



/(*) = /„ + £ 15^ 



i=l 



180) 



Thus, the multi-center metric admits multi-NUT singularities so that the jth NUT singular- 
ity is characterized by the charge rrij and is located at Sj . However, we can remove these 
singularities provided all NUT charges are equal, 



mi = m 2 



mi 



9 



In this case, the cyclic variable x 4 is periodic with the range, 

9 



< x 4 < 4tt 



N 



(181) 



;i82) 



We can now investigate the symmetries associated with the projection of the particle's 
motion onto the curved 3-manifold described by the multi-center metric tensor, 

N 



9jk{x) 



i=i 



;i83) 
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The projected Hamiltonian is given by 

1 _ 2 

T-i = -IT 2 + f(x)W(x) with W(x) = U(x) + jf(x) + £f~\x) - £ . (184) 



Let us first note that for multi-center metric (183), it is straightforward to deduce from 



Theorem III. 2 the metric condition, 



E(n ■ x)di — (n ■ 3i)x 
— — : = U 



(185) 



which can not hold for more than two centers. Thus, we state the following theorem 
Ngome 08/20091 : 



Theorem III. 4. In the curved 3-manifold carrying the N '-center metric, 



N 



9jk{x) 



/o+E 



rri; 



jk , 



no symmetry of the Kepler-type occurs for N > 2. 

For simplicity, from now on we limit ourselves to a discussion of the two-center metrics, 

(186) 



f(x) = f + ^r 1 ^ + , \x±a\^o, 



\x — a\ \x + a 

which are relevant for diatomic molecule systems and which possess some interesting sym- 
metry properties. These metrics include, as special regular cases, those listed in Table [Tj 



/o 


N 


Type of Metric 





1 


(mi or m,2 = 0) Flat space 


1 


1 


(mi or m 2 = 0) TAUB-NUT 





2 


Eguchi-Hanson 


1 


2 


Double TAUB-NUT . 



TABLE I: Examples of two-center metrics. 



Let us now apply the van Holten algorithm ( 136 ) to derive the symmetry of the two-center 
metric. 
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1) First, finding conserved quantities linear in the covariant momentum require to solve 
the reduced series of constraints, 



C m d m (f(x) W(x)) = o(0) 
d n C = qF nm C m o(l) 
V i C l +V l C i = 0. o(2) 



(187) 



From Theorem III. 1 , we deduce the rank-1 Killing tensor satisfying the second-order 



constraint of (187), 



Ci — 9im{%) £ Ik X i a — \ 



188) 



The tensor (188) generates rotational symmetry around the axis through the two centers. 



Next, injecting both (188) and the magnetic field of the two centers, 

x + a 



-* X — Qj 

B = mi — — — + m 2 
\x — ar 



into the first-order equation of (|187|) yield 



q mi 



\x — a\ 



+ m 2 



\x + a\ 



x + a 
\x + a\ 



a 
a 



(189) 



(190) 



Finally we obtain, as conserved quantity, the projection of the angular momentum onto the 
axis of the two centers, 

x + a 



J a = C a - q ( mi — % + m 2 



\x — a 



x + a 



with L a = [x x II • 



(191) 



which is consistent with the axial symmetry of the two-center metric. 



2) Now we study quadratic conserved quantities, Q = C + C l Hi + - C' lJ iTIL- . Putting 
Cijk = Cijki = • • • = , leaves us with, 

' C m d m (f(x) W(x)) = o(0) 

d n C = qF nm C m + C n m d m (f(x)W(x)) o(l) 

Vid + VtQ = q (F im C x m + F lm Cr) o(2) 

^ ViQj + VjCii + V^j = . o(3) 
We consider the reducible rank-2 Killing tensor, 



(192) 



a 



2 2 

~2 9im(x) 9jn(x) e m [k e" pg a 1 a? x k x q + — g a (x) 9j m {x) a 1 a 



l „m 



a- 



(193) 
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which is a symmetrized product of Killing- Yano tensors. C, = gi m (x) e 



ira V 
m 



lk — x generates 
a 



rotations around the axis of the two centers and Cj = gj m (x) — generates spatial trans- 
lation along the axis of the two centers. Injecting (193) into the second-order constraint of 



(192) yields 



2q 

a? 



d = 1 g im e m jk a 3 x k ai ( ;//., 



x l — a 1 



\x — a 



+ m 2 



x l + a 1 
x + a\ 



For vanishing effective potential, we solve the first-order constraint with 

2 



c 



q 



(x l — a 1 ) (x l + a 1 ) 

mi A 7 ^ — ^— a\ + m 2 \^ — ^- a x 



\x — a\ 



\x + a 



(194) 



(195) 



so that we obtain the square of the projection of the angular momentum onto the axis of 
the two centers, plus a squared component along the axis of the two centers of the covariant 
momentum, 

Q = JI + K- (196) 

As expected, this conserved quantity is not really a new constant of the motion 
[Gibbons 19 871 |Valent 09/2003[ |Valent 07/2004[ [Duval 05/2005] . 

• Now we turn to the Kepler-type dynamical symmetry. Let us first check if a rank-2 
Killing tensor associated with Runge-Lenz type conserved quantity does exist. To this end 
we apply Theorem III.2| to the two-center metric, 

mi rn 2 



We obtain 



9jk(x) = f(x)5 jk , f(x 
n x (x x V f{x) 



\x — a\ 



+ 



\x + a 



m 2 



mi 



x + a 



\x — a\ 



(x x a) x n = 



according to Theorem III. 2 Consequently we get 
m 2 Tti\ 



\x + a\ 







or 



x = ka , k = const . 



\x — a\ 



(197) 
(198) 

(199) 



The right condition in (199) restricted to motions parallel to a and therefore implies no 
interesting case. 



Considering the first case given by (199), we assume that both charges are positive mi > 
, m 2 > , and we write a = (aj., a 2 , a 3 ) . Thus, the left Equation in ( |199 ) becomes 



(x 



a iP) +(y-a 2 p) +(z-a 3 p) 
with p 



(P 2 ~ 1) 



2/3 . 2/3 

mi + m 2 



(200) 



2/3 



m 2 — m l 



2/3 
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We recognize here the equation of a 2-sphere of center dp and radius R = a \/ p 2 — 1 , noted 
as S 2 . The latter shows that for two-center metric, a Kepler-type dynamical symmetry is 
only possible for motion confined onto the sphere, S 2 . 

Before searching for the exact form of the associated Runge-Lenz conserved quantity, let 
us first check that the motions can be consistently confined onto this 2-sphere, S 2 . 

To this end, we assume that the initial velocity is tangent to S 2 and, using the equations 
of motion, we verify that at time (t + 6t) the velocity remains tangent to S 2 . Thus we 
write 

v(t + St) = v + 5t v with v = v(t ) tangent to S 2 . (201) 



The equations of motion in the effective scalar potential chooses as (205) have the shape 

,2 



U = qv x B - V {f{x)W{x)) 



V 



mi 



\x — a\ 



+ 



m 2 



\x + a\ 



x . 



(202) 



Injecting the expressions of the magnetic field of the two-center (189) and the effective 



potential ( |205[ ), we obtain 

m 1 



m 2 



\x — a\ 



\x + a\ 



rrii 



\x — a\ 



+ 



m 2 



\x + a 







x = f(x)vo. 



Thus v(to + 5t) in (201) becomes 



nil 



+ 



m 2 



q 



nil 



+ 



m 2 



\x — a\° \x + a| 3 J \* \\x — a\ \x + a\ 

where vq and xq are tangent to the 2-sphere S 2 . Hence, the velocity remains tangent to S 2 
along the motion. 

Having shown the consistency of motions on the 2 -sphere S 2 , we can state the following 
theorem Ngome 08/2009] : 



Theorem III. 5. In the curved 3-manifold carrying the 2-center metric, 

nil rn2 



9jk{x) 



fa 



+ 



\x — a 



\x + a 



(203) 



a scalar Runge-Lenz-type conserved quantity does exist only for a particle moving along the 



axis of the two centers or for motions confined on the two-sphere of radius R = a y p 2 — 1 
centered at a p (mi, m 2 > 0) . In the Eguchi-Hanson case (mi = m 2 ) , the 2-sphere is 
replaced by the median plane of the two centers. 
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Our method here is particular since instead searching for Kepler-type dynamical sym- 
metry directly, we already look at the conditions of its existence. Knowing now, for the 
two-center metric, that only motions confined on S 2 allow a Runge-Lenz-tpe conserved 



quantity, we can solve the second-order constraint of (|192|) using (|200|). Thus, we obtain 

Q 



C; 



jk 



mi 
\x — a\ 



+ 



m 2 
\x + a\ 



(204) 



where the only component of n is along the axis aj a . The final step consists to solve the 



first-order constraint of (192). Indeed, following (161) the clue is to choose 

2 



f(x)W(x) 



q 



mi 



\x — a\ 



+ 



m 2 



\x + a\ 







Till 



\x — a\ 



+ 



m 2 



\x + a\ 



+ 7- 



(205) 



with (3, 7 G R . Let us precise that this potential satisfies the consistency condition given 



by the zeroth-order constraint of (192). Moreover, the leading coefficient of the effective 
potential cancels the obstruction due to the magnetic field of the two centers, and the 



remaining part on the right-hand side of (|205|) leads to 

x + a 



C 



> / dC a 
P mi — — — + m 2 
\x — a\ 



x + a 



a 
a 



(206) 



Collecting our results (144), (204) and (206) provide us with the scalar, 



K n 



Tlx J 



a , „ „ 
- + -(C a -J a 
a q 



(207) 



which represents, in the case of two-center metrics ( 186 ), a conserved Runge-Lenz-type scalar 



for particle moving on the 2-sphere of center positioned at a p and radius R = a \/ p 2 — 1 



combined with the effective potential (205). 
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D. Killing- Stackel Tensors on extended manifolds 



Having discussed, in the two previous sections, the conditions on Killing tensors which 
are related to the existence of constants of motion on the dimensionally reduced curved 
manifold. We can observe that the Killing tensor generating the Runge-Lenz-type quantity 
preserved by the geodesic motion can be lifted to an extended manifold. 

Let us study this lifting problem in detail, by considering the geodesic motion of a particle 



before dimensional reduction ( 129 ). The particle evolves on the extended 4- manifold carrying 
the metric g^x) with v — 1, - • - , 4 . A rank-2 Killing-Stackel tensor on this curved 4- 
manifold is a symmetric tensor, Cn U , which satisfies 







A, /i, v = 1 



(208) 



For the Killing-Stackel tensor generating the Runge-Lenz-type conserved quantity, the 
degree-2 polynomial function in the canonical momenta p^ associated with the local co- 
ordinates x^ , 

K= l -C^p, Pu (/!, v = !,■■■ ,4), (209) 

is preserved along the geodesies. Then, the lifted Killing-Stackel tensor on the 4-manifold, 
which directly yields the Runge-Lenz-type conserved quantity is written as 




i, i = 1,2,3. 



(210) 



The tensor Cjj is, therefore, a rank-2 Killing tensor on the dimensionally reduced curved 3- 
manifold carrying the metric gij{x) = f(x) Sij , which generates a Runge-Lenz-type quantity 
conserved along the projection of the geodesic motion onto the curved 3-manifold. The off- 
and the diagonal contravariant components read 

1 



C 



M 



c 



q 



-C l - C\ A K , C 



1 1 



(2/q 2 )C - (2/q)C k A k + C jk A j A h 



(211) 



The term A k represents the component of the vector potential of the magnetic field. In the 



case of the generalized TAUB-NUT metrics, the terms C and are the results (474) and 



(468) of the first- and the second-order constraints of (156), respectively. In the case of the 



two-center metrics, C and are given by the results (206) and (204), respectively. 
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As an illustration, let us consider a particle in the gravitational potential, V(r) 



m G 



described by the Lorentz metric |Bargmann Duval 05/1984, IBalachandran 1986|, 
IDuval 1991 j . 

dS 2 = dx 2 + 2dx 4 dx 5 -2V{r) (dx 5 ) 2 . (212) 

The variable x 5 = t is the non-relativistic time and x 4 the vertical coordinate. Rotations, 
time translations and "vertical" [on the fourth direction] translations generate as conserved 
quantities the angular momentum L , the energy and the fixed mass m , respectively. The 
Runge-Lenz-type conserved quantity, along null geodesies of the 5-manifold described by the 



metric (212), 



K =\ C ab PaPb with a, b, c = 1, 



is derived from the trace-free rank-2 Killing-Stackel tensor [Duval 19 91 J. 



C 



ab 



g ab - rj ab 



with f) = r] ab g. 



ab 



For some n G R 3 , the nonvanishing contravariant components of r\ are given by 



n l x J + n J x % — fj 8 %3 and r] 



15 



n 
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V = riiX , 



(213) 



(214) 



(215) 



where we recognize, in the left hand side of (215), the generator of Kepler-type symmetry 
in the dimensionally reduced 3-manifold. 



A calculation of each matrix element of the Killing tensor (214) leads to C whose only 
nonvanishing components are, 

C ij = 2 fj S ij - n l x j - n j x i and C 44 



2f)V{r) . 

Consequently, the associated Runge-Lenz-type conserved quantity reads as 



K ■ n = [p x L + m V(r) x)-n. 



(216) 



(217) 



Note, in the previous expression, that the mass " m " is preserved by the "vertical" reduction. 
In the original Kepler case, we thus deduce on the dimensionally reduced flat 3-manifold 
that the symmetric tensor 

C ij = 2 5 ij n k x k - n L x j - n j x i (218) 

is a Killing-Stackel tensor generating the Runge-Lenz-type conserved quantity along the pro- 
jection of the null geodesic of the 5-manifold onto the 3-manifold carrying the flat Euclidean 
metric. 



53 



IV. NON-ABELIAN GAUGE FIELDS AND THE BERRY PHASE 



Conserved quantities of an is ospin- carrying particle in non-Abelian monopole-like fields 
are investigated. In the effective non-Abelian field for nuclear motion, obtained through the 
Berry phase in a diatomic molecule, due to Wilczek et ai, an unusual conserved charge and 
angular momentum are constructed. 



A. The Wu-Yang monopole 

In 1968, T.T. Wu and C.N. Yang found the first "monopole-like" classical solution of 
the Yang-Mills field equations |Wu Yang 1968 . Such a solution can also be viewed as an 



extension of the Abelian Dirac monopole solution when usual electrodynamics, with £7(1) 
symmetry, is considered as a part of a larger theory. The generator of the electromagnetism 
£7(1) subgroup should be embedded into the non-Abelian SU(2) gauge group. 

In order to investigate the Wu-Yang monopole solution, we consider here a Yang-Mills 
theory described by the Lagrangian density, 

C = -\F« u F a ^, //,!/ = 0,1, 2, 3. (219) 

where represents the antisymmetric Yang-Mills field strength tensor taking values in 
the Lie algebra of the gauge group SU (2) , 

F; v = d,A a u - d v Al - e e ahc A\Al , a = 1, 2, 3 . (220) 

Here e and the antisymmetric tensor e a bc denote the gauge coupling constant and the 
structure constant of the gauge group, respectively. As expected, 

A, = Alr\ (221) 

takes value in the su(2) Lie algebra. The Hermitian and traceless infinitesimal generators 
of the SU (2) gauge group verify the commutation relation, 

[r a ,r b ] = %e ahc r\ r a = l -a\ (222) 

where the a a are Pauli matrices. Our approach now is to find the equations of the motion 



using a variational principle with the Yang-Mills Lagrangian density (219). We obtain the 
classical source-free Yang-Mills (YM) equations 

d p F daP - ee adc A% F aafi = , (223) 
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which can be written in a more compact way as 



DpF iaP = 0. 



(224) 



Searching for solutions of the equations of the motion (224) in the "temporal" gauge 



A a = 



(225) 



only time-independent gauge field and local gauge invariance are permitted. We can now 
posit the spherically symmetric Wu-Yang ANSATZ |Wu Yang 1 968 1, 



-Aj g Eiaj X' 



;i-$(r)) 



(226) 



r 2 = XiX 1 with i, j = 1,2, 3 



where a, and g represent the color index, the space indices and the quantized Wu- 

Yang monopole charge, respectively. Also remark that $ is a radial real function which is to 
be determined. As expected, a direct calculation implies that all time-dependent components 
of the field strength vanish, 

F o ; = 0, // = 0,1,... ,4, (227) 
whereas the spatial components of the 2-form curvature reads 



F$ = ge ijk {2fi 



«*(!-*) 



k^a 



— egx x 



1 - $ 



g d (\ — <3>\ , i l j i 

H r~ I I ] £jal XX &ial X X f . 

r ar \ r z 1 



(228) 



The right hand side bracket of (228) can be rewritten by using the relation, 



(e jal x i x l - e ial x j x l )r a = e ijk (r 2 5 ak - x a x k )r a 



(229) 



so that the 'Wu-Yang" field strength reduces to 



8 ak d$ x a x k (d& $ 2 , N $ 2 -eg 

— — + —3- — - eg— + 2(eg-l)-+ y - 

r ar r 6 \ar r r r 



(230) 



Injecting the relations (226) and (230) into the Yang-Mills field equations (224) provide us 



with the non-linear Wu-Yang equation, 
„ d 2 $ 



dr 2 



eg(l - eg + eg$) ($ - l) ( $ + 



2 -eg 

eg 



0. 



(231) 
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Note that due to the non-linear nature of the YM equations, to search analytical solutions 



of (231) is an unconquerable task. 



Hence, we investigate numerically the non-linear equation (231) viewed as a dynamical 



system. To this end, without loss of generality, we reduce ourselves to the case where eg = 1 
2 . Thus, the equation (231) takes the simple form, 

„ d 2 $ 



dr 2 



$($-l)($ + l) =0. 



(232) 



We first posit the variable change 



r = exp(r) , r e R, r E IT 1 



(233) 



where r i s vie wed as an evolution parameter. Next, we multiply the resulting r -dependent 

J I d$(r) , 
equation (232) by — - — so that we get 



dr 



d 1 fd$(T 



dr \ 2 \ dr 



4 



(* 2 (r)-iy 



dr 



(234) 



Then, the equation (234) can be interpreted as the equation of motion of a unit mass particle 
with non-conserved Hamiltonian 3 , 



H=^$ 2 + V with V = -- ($ 2 - lV 
2 4 v ' 



(235) 



Here V represents an Higgs potential in which the particle evolves. Let us remark that the 
kinetic frictional force, 

F + = $, (236) 

exerted on the particle has a positive friction coefficient and makes the energy to grow, since 

dE 



dr 



> 



(237) 



The system receives energy from the exterior so that the Rayleigh function, 1Z , is negative 
and reads 



n 



(238) 



o 



III A 



2 See the formula ( 94 ) in section 

3 Here the dot means derivative w.r.t. the evolution parameter 



d_ 
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When positing the conjugate momentum as $ = $ , we can construct the canonical phase- 
space (<&, \fM in which we define an extension of the equations of the motion of our non- 



conservative system (234) as 




(239) 



* + $($ 2 - l) . 



Then, we can now describe the curve solutions of this Hamiltonian system (239 ) Protogenov[ 
IBreitenlohnerj . To this end, we draw in the phase-plane (<&, \&) , the vector field ($, , 
representing the velocity of each phase-point. Hence, the orbit solutions of the dynamical 
system lie on curves tangent to the velocity vector field. However, we restrict our inves- 
tigation to finite orbits which are the only solutions physically consistent. We first search 
for critical points, which can be considered as orbits degenerated to a point, by solving the 
constraints 

' V = 0, 

(240) 

^ + $($ 2 - 1) =0. 

A simple algebra leads to the three critical points ($, = { (0, 0) ; ( ± 1, 0) } , to which 
we characterize the equilibrium by analysing the eigenvalues and the eigenvectors of the 
stability matrix, 

d 2 n \ 

I U 1 

(241) 



/ 



A 



o 2 h 



d^ 2 
d 2 H d^R 
+ ) 



3$ 2 - 1 



/ 



9$ a& 9$ 

For the fixed point ($, = (0, 0) , the eigenvalues of the stability matrix read, 

l/2±i>/3/2. 



A 



± 



As the complex conjugated eigenvalues \± have both a positive real part, then the orbits are 
spiraling out with respect to the focus (0, 0) . Hence, the critical point (0, 0) is unstable 
and is considered to be a negative attractor. 

• For the fixed points ($, \fn = ( ± 1, 0J , the eigenvalues of the stability matrix A read 



A, 



2 and A_ 



-1. 
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The fixed points are saddle points with the stable direction given by A_ and the unstable 
direction given by A + . 



• We also consider the two bounded solutions, noted $=|= , represented by the curves joining 
the negative attractor (0, 0) to the saddle points (=Fl, 0). The corresponding curves solution 
are represented in the phase portrait below by $_ in pink line and $ + in cyan. 




FIG. 5: Phase portrait of the Hamiltonian system (239). 



Let us now discuss the bounded solutions of the non-linear differential equation (232). 
For r e R + , we have the five orbits solution, 



$(r) = {-1, 1, 0, $+} 



(242) 



that we introduce in the Wu-Yang forms (226) and (230) to obtain the shape of the gauge 
field. 



The degenerate orbit solution $ = — 1 leads to the pure gauge field 

A? = 2ge iaj - since = F? k = . 
The choice $ = 1 corresponds to the null gauge potential, 

A$ = with F£ = i^ = 0. 
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(243) 



(244) 



It is worth noting that the two previous gauge potentials, with both vanishing field strength, 
can be transformed the one by the other using a suitable gauge transformation. They are 
therefore gauge equivalent. 

• We consider the two curves solution $_ and $ + , admitting the asymptotic limits, 

$^ when r > 1 f lim $± = ±1 

$± = { ' with I ™ (245) 

$° ± when r < 1 [ lim $^ = . 

Taking into account the behavior of the curves solution in the neighborhood of the origin, 
we can neglect the cubic term of the non-linear equation (232). Thus, we deduce that $![ 
satisfy the differential equation, 

r 2 d ±} J +(f>l(r) = 0, (246) 
which provides us with the non-analytic solution, 

$o = ± ay /rcos | —In — | , with (a, r < l) G R* x R + * . (247) 
\ 2 r / 



In the case where r> 1, the equation (232) reduces to the simple form 



#^(r) const /I 



<ir 2 r 3 

so, we derive the behavior of the radial functions Q± at infinity as 

$°?(r) = ±1 =f ^ + O f -^J , where 7 > . (249) 

After the complete description of the asymptotic behavior of the functions <&± , our business 
now is to fill the gap between these two limit cases. Here we investigate the "solution" K + 
but the case K_ is not more complicated. Thus, we integrate numerically the non-linear 



equation (232) from a point ro , located in the neighborhood of the origin so that $ + is 
approximated by 3>5_ , till a sufficiently great value of r so that $ + can be approximated by 
. Following the usual procedure, let us begin the numerical integration by adding some 
analytic correction terms, , into the early expression (247) of 3>9 . Thus, we express the 



numerical lower bound as 



$°(r) = iav^cos ( ^ln- ) +Va i («, r ) r\ (250) 
1 2 r 0y / ^ 
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where the coefficients a« depend on the values of ro <C 1 and the fixed parameter a . For 
the fixed values of the integration parameters, 

r = 0.007873997658 , a = -0.8873554901 , (251) 

and with the initial conditions of the numerical integration given by 

$° ± (r ), ^i(ro) 



we obtain, for r G [ro, 14] , the curve solution $+(r) of the non-linear equation (232). 



For e — g — 1, we draw the field strength intensity of the usual SU(2) Dirac monopole 



(in circling dashes) together with the intensity of the field strength solution (230), namely 
B±(r) — — (l — $^_) , carrying the branch $ + (in heavy line). 

It is worth mentioning that, in the asymptotic limits, the two curves coincide. The latter 
is due to the fact that when the radius tends to zero, B±(r) corresponds to the length of 
a Dirac monopole field; and when the radius tends to infinity, we get the null length of 
the vacuum. In the case of intermediary values of the radius, the two curves become quite 
different. This is a consequence of the radial function $ + which takes non-quantized but 
continuous values between zero and one. (See Figure [6]). 

• Let us now investigate the last bounded solution, Q± = . The comparison of the 
intensities of field strength made above provides us with strong assumption on the nature 
of the solution when <&± = 0. Let us analyze further by injecting the solution $(r) = 



into the Wu-Yang ansatz (226). In that event, we obtain the gauge field 



A* WY = ge m ~, (252) 



which implies that the field strength (230) reduces to that of a Wu-Yang monopole 



The energy density, 



Ftr = 9^~^- ( 2 53) 



m = \F% W¥ F aijWY = £;, (254) 



is singular at the origin r = , so the Wu-Yang monopole possesses an infinite magnetic 
field energy. 

Moreover we can prove, as suggested by the figure in the bottom of (|6]), that the non- 
Abelian Wu-Yang monopole can be viewed as an imbedded Dirac monopole. Indeed, we 
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FIG. 6: In the top side we plot the curve solution <£ + ; and in the bottom side we compare the 
field strength intensities of B + with the imbedded Dirac monopole. 



consider the trivial imbedding of Dirac monopole field A^ u ^ into £577(2) , 

a d u(i) _^ a d sum = A« D r a , (r a = a J 2) , (255) 

where the Abelian gauge potential reads 

for a =1,2, 

K° = { ( 256 ) 
Al D = ±g(l T cos9)d^. 
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Thus, a short algebra yields the imbedded gauge potential, 



A 



A DSU(2) 



DSU(2) 




(temporal gauge) 



aDSU(2) , (l^COS^) . 

A- = ±g ^ — - — (— sm0, cos0, 0) r 3 



(257) 



rsm9 

In order to avoid the Dirac string singularity, we apply the singular "hedgehog" gauge 
tr ansf ormat ion , 

A DSU(2) _^ u f A DSU(2) + = A WY ^ (25g) 

which rotates the unit vector on the sphere S 2 to the third axis in isospace. Thus, U is 
charaterized by the unitary matrix, 

6 9 



U{6, 



cos 



sin ^ exp (inc, 



— sin - exp [— m<; 







cos 



/ 



n G IN*. 



(259) 



Applying (|258|) with (|259|), it is straightforward to derive the gauge equivalent potential, 

(temporal gauge) , 



AaWY 
^0 



.4 



1WY 



71 



cos 9 sin 9 cos (nd)) did) H — sin (n0) 9^ . 

e e 



2 



n 



cos 9 sin sin (n</>) did) cos (nd)) di9 , 



(260) 



A 



3WY 



71 



sin 2 #a</>. 



When taking into account the Dirac quantization relation 4 , eg = n = 1 , and the following 
algebraic relations, 



d,,9 



x 



1 

, , zcos(p, zsmif,, 
r z V cos ( 



and dj( 



1 



r 2 sin 2 9 



-y, x, 0) 



(261) 



AaWY 
^0 



the gauge potential p60j) transforms into the cartesian form as 

(temporal gauge) , 

A\ 
A? 



IWY _ 

i 


£(<>. 


z, -y) 


2WY __ 

i 




z, 0, x 


3WY 

i 




-x, 0) 



(262) 



See the formula (94) in section III A 
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By compacting (262), we hence recover the exact expression [see (252) and (253)] of the 
non-Abelian Wu-Yang gauge potential with a "hedgehog" magnetic field. It is now clear 
that the non-Abelian Wu-Yang monopole field can be obtained by imbedding the Abelian 
Dirac monopole field into an 577(2) gauge theory. 

Let us now inquire about conserved quantities. To do this, we first Identify the su(2) Lie 
algebra of the non-Abelian generator, r a , with R 3 . Thus, we make the replacement, 



r a — y X a , with the internal index a = 1,2,3 



(263) 



Here the non-Abelian variable X a represents the isospin vector which satisfies the Poisson- 
bracket algebra, 

{X a ,X b } = -e abc X c . (264) 

Hence, we consider an isospin-carrying particle [Balachandran 1977| IDuval 19781 
IDuval 19801 IDuval 1982] moving in a Wu-Yang monopole field |Schechter[ IBoulware 19761 
IStern 19771 |Wipf 1986 , augmented by a scalar potential |Schonfeldl IFeher 1984] . described 
by the gauge covariant Hamiltonian, 

?2 



7T- 



H = — + V(x,X a ), n=p- eX a A 



a 7a WY 



(265) 



We define the covariant Poisson-brackets as 



dg df 



{f>9} = DifZ--^D j9 + erF?™- 



df dg _ df dg 

£abc r* — ~ — ' 



8lTj dlik 



3X a dX b ' 



(266) 



where Dj is the covariant derivative, 

Djf = djf - ee ahc X a A) 



|/< WY df 

dX c 



(267) 



Thus, the commutator of the covariant derivatives is recorded as 

d 



[A, Dj 



—€abc£ a FAA 



b WY 



(268) 



13 dX c ' 

Following van Holten's recipe |van Holten 2007] . conserved quantities Q[x, X, 7f) can conve- 
niently be sought for in the form of an expansion into powers of the covariant momentum, 



Q(x,X, ?F) = C(X,X) + Ci(x,X)TTi + —Ci j (x,X)TT i 7T j + 



(269) 
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Requiring Q to Poisson-commute with the Hamiltonian, {<2,"H} = , provides us with the 
set of constraints to be satisfied, 
BC BV 

W + e ak l^-=0, o(0) 

BC BV 

D t C = eX^F^C, + CyDjV + e^-^—, o(l) 

BC ■ BV 

D i C j + Djd = eZ a {F? k WY C k3 + Ff k WY C kl ) + C ijk D k V + e abc X a ^ — , o(2) 

DiCjk + DjCki + DkCij = eX a {F°i WY Cijk + F°i WY Ciki + F kl WY Cuj) 

.dC ijk BV 



+CijkiDiV + e abc X a - 



BX b BI C 



o(3) 



(270) 



To start, we search for zeroth-order conserved quantity. Thus C{ = C, 



'j 



, so that 



the series of constraints ( 270 ) reduces to 



BC BV 
DiC = . 



(271) 



o(l) 



The zeroth-order equation of (271 ) is identically satisfied for rotationally invariant potentials 



with respect to x and X . Applying for the consistency condition, 

BC 



[D u D J ]C = -e abc X a F^ WY — 







we get the shape of the derivative of C along the isospin variable, 

BC 

f(r,X)x c + h(r,X)X c . 



BX C 



(272) 



(273) 



Injecting (273) into the first-order constraint of (271) implies, for an arbitrary function 
/i(r, Z) , that 

B % C = f(r,X) (xi - , [with eg = l]. (274) 

From the commutation rule, 

[B i ,B j ]C = I i 



r \ar r 



(275) 



we derive the exact form of the function f(r,X) 



f(r,X) 



(3 = const G E . 



(276) 
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Taking into account the result (276), the equations (274) and (273) lead to the system of 
equtions, 



diC = p 



dC 

~dT a 



(277) 



which is solved uniquely by the covariantly constant charge, 



Q(x,X)=/3 



x ■ X 



7 g(Z), 7 = const GR. 



(278) 



This charge can be viewed as a linear combination of two quantities separately conserved 
along the particle's motion since /3 and 7 are arbitrary real numbers. Note that the first 



term on the right-hand side of (278), namely, 



x ■ X 



(279) 



can be seen as a conserved electric charge; and its conservation admits a nice interpretation 
in terms of fiber bundles Horvathy 12/1984, Horvathy 06/198"5] . The SU(2) gauge field is 
a connection form defined on a bundle over the 3-dimensional space so that for Wu-Yang 
monopole, the su(2) connection living on the (trivial) bundle reduces to the U(l) Dirac 
monopole bundle. This is the reason why the electric charge is conserved in the Wu-Yang 
case: the latter is, as already seen, an imbedded Abelian Dirac monopole. 

Next, we study conserved quantities which are linear in 7Tj, CV,- = Cjj& = ■ ■ ■ = 0. We 
therefore have to solve the constraints, 



( dC dV 

C l D i V + e abc l a — h —=0, 



ra rpa WY , 



o(0) 



f)n. pj\r 

D i C = el"F^C j + e abc l"p b ^ c ,o(l) 



(280) 



k DiCj + Djd = . 



o(2) 



When the potential is invariant with respect to joint rotation of x and X, inserting the Killing 
vector generating the spatial rotations, 



C = n x x . 



into the series of constraints (389) yields 



C = -Q 



n ■ x 



(281) 



(282) 
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Collecting the two previous results provides us with the conserved angular momentum, 

x 



J = X X 7T — Qc 



Let us now turn to quadratic conserved quantities, Cj 



ijk 



a 



ijkl 



(283) 
0. The set of 



constraints (270) reduces to 



C^V + e abc ^— b — = 0, 



DiC = eZ a F* WY C 3 + C %3 D 3 V + e ahc T 



da dv 



dX h dX c '' 

DiCj + DjCi = eX a {F° k WY Ckj + Fj k wr Cki) + t a b c T 
DiCjk + DjCki + DkCij = . 



o(0) 
o(l) 



(284) 



o(3) 



We observe that the rank-2 Killing tensor generating the Kepler-type dynamical symmetry 
has the property, 

Cij = 25ij n ■ x — (riiXj + n 3 Xi) . (285) 



Inserting (285) into (284), from the 2nd-order equation we find, therefore, 

C=Qo— • 



(286) 



The first-order equation requires in turn 

DiC = \ Ql (n ■ x) ^ - Q 2 -| \ + 2 n ■ x D { V - n { x ■ DV - x^n ■ DV . (287) 



Restricting ourselves to potentials falling off at infinity and invariant by rotations with 
respect to x and X , we hence obtain 



n ■ x 



V = + - + (3 and C = a 

2r r r 



(288) 



where a and /3 are arbitrary constants. Note that the inverse-square term in the previous 
potential is fixed by the requirement of canceling the bracketed term on the right-hand side 



of (287). Collecting our results yields, 



x 



K = 7? x J + a- , 



(289) 



which is indeed a conserved Runge-Lenz vector for an isospin-carrying particle in the Wu- 



Yang monopole field, combined with the potential (288) |Horvathy 1991| . Let us emphasize 
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that the fine-tuned inverse-square term is necessary to overcome the obstruction in solving 
the constraint equation Ngome 02/2009 ; without it, no Runge-Lenz vector would exist. 



The expression (289) is actually not surprising since the Wu-Yang monopole field, as we 
proved it, corresponds to an imbedded Dirac monopole field. 

The importance of the conserved quantities J and K is understood by noting that they 
determine the trajectory : multiplying the conserved angular momentum by the position, 
x, yields 

J-- = -Qo, (290) 

r 

so that the particle moves, as always in the presence of a monopole, on the surface of a cone 
of half opening angle, 

ip = arccos(|Q |/^) where J = \J\ ■ (291) 
On the other hand, the projection of the position onto the vector N, given by, 

N = K + {a/Q )J, N-x = J 2 -Q 2 = const, (292) 

implying that the trajectory lies on a plane perpendicular to N. The motion is, there- 
fore, a conic section. Careful analysis would show that the trajectory is an ellipse, a 
parabola, or a hyperbola, depending on the energy, being smaller, equal or larger than 
(3 jFeher 19841 IFeher 19851 IFeher 19861 IFeher 19871 ICordani 1990] . In particular, for suffi- 
ciently low energies, the motion remains bounded. 

The conserved vectors J and K satisfy, furthermore, the commutation relations 

{Ji, Jj] = e ijk J k , {Ji,Kj} = e ijk K k , {K u Kj] = 2(/3 - H)e ijk J k , (293) 

with the following Casimir relations, 

J-K = -aQ , K 2 = 2{U - P){J 2 - Ql) + a 2 . (294) 

Normalizing K by [2{f3 — %)] 1 , we get, therefore an SO(3)/E(3)/SO(3, 1) dynamical 
symmetry, depending on the energy being smaller /equal/larger than /3 5 . 

We remark that although our investigations have been purely classical, there would be 
no difficulty to extend them to a quantum particle. In the self-dual Wu-Yang case, the 

5 The dynamical symmetry actually extends to an isospin-dependent representation of SO (4, 2) 
|Horvathy 1991 . 
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SO(4)/SO(3, 1) dynamical symmetry allows, in particular, to derive the bound-state spec- 



trum and the Scattering-matrix group theoretically, using the algebraic relations (293) and 



(294) jFeher 19841 [Cordani 19881 IFeher 1988 1 IFeher 1 989J. 
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B. The Berry phase - general theory 



The Berry phase arises for systems which can be conveniently described in terms of two 
sets of degrees of freedom |Berry 1984 . The one is "fast" moving with large differences 



between excitation levels, and the other is "slow" with small associated energy differences. 
This decomposition is extensively used in molecular physics through the adiabatic or Born- 
Oppenheimer approximation. In a molecule, for instance, the electronic motion is described 
by the "fast" variables f , and the nuclear motion by the "slow" degrees of freedom R . 

We first deal with the "fast" degrees of freedom, keeping the "slow" as approximately 
fixed. In that event, we simply solve the stationary Schrodinger equation for the "fast" 
variables, with the "slow" variables appearing parametrically, 

n{R)^ m {r, R) = S m {R)y m {r, R) . 

Next, we complete the analysis by allowing slow variations in time for the previously fixed 
variables. The (adiabatic) assumption is that the slowly varying degrees of freedom R do 
not change quickly enough to induce transitions from one S n level to another. Thus, the 
system starting in an initial eigenstate remains in this state in response to the slow change 
of the variables R appearing parametrically. 

As a consequence of this effective dynamics, an external vector potential called the 
Berry connection is induced. It has been argued that the "feed-back" coming from the Berry 



phase modifies the (semi) classical dynamics |Chang 1995 INiu] . The associated magnetic 



like field is the Berry curvature, and the line integral of the connection is the "celebrated" 



Berry phase |Berry 1984 , ISimonj . See also |Aitchison] . 



For better understanding, let us study deeper the way to obtain the Berry gauge poten- 
tials. To this end, we separate the following in Abelian and non-Abelian cases. 



The Abelian gauge potential : non-degenerate states 



Following Berry's original paper [Berry 1984| , let us point that an U(l) gauge field 
may appear when a single non-degenerate level is subject to adiabatically varying external 
parameter. To this end, we consider a physical system described by a Hamiltonian which 
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depends on time through the vector R(t) , 

H = H(R) , R = R(t) . (295) 

Here R(t) denotes a set of m classical parameters, 

R(t) = (R 1 (t),R 2 (t),--- ,R m (t)), 

slowly varying along a closed path C in the parameter space, since the system is assumed to 
evolve adiabatically. We introduce an instantaneous orthonormal basis constructed with the 
eigenstates of H(R(t)) at each value of the parameter R. The eigenvalue equation reads 

H(R(t))\m 1 R(t)) = S m (R(t))\m 1 R(t)). (296) 



The basis eigenfunctions \m,R(t)) are not completely determined by (296). Indeed, without 
loss of generality, while normalizing the eigenfunctions, 

(m , R(t)\m , R(t)) = 1, 

this implies that eigenfunctions are unique up to multiplication by a phase factor. More- 
over, for a slowly varying Hamiltonian, the quantum adiabatic theorem states that a system 
initially prepared to be in one of its eigenstate \m , R(0)) , at t = , remains in this instan- 
taneous eigenstate along the cyclic process C. Consequently, the quantum state at time t 
can be written as, 

\^ m (t)) = exp (ia m (t)) \m , R(t)) , (297) 



where the exponential term in (297) is the only degree of freedom we can have in the quantum 



state. Substituting the expression (297) into the Schrodinger equation 



ih^\* m {t)) = H(R(t))\*m(t)) , (298) 
and projecting both sides of the equation onto (n , R(t) \ , yields the equation for the phase, 

a n = i I \m , R\V R \m , R) • dR- \ [ £ n {R{t'))dt' . (299) 
Jc " n Jo 

Thus, in addition to the usual dynamical phase, 



ijf E n {R{t'))dt' , (300) 
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the quantum state acquires an additional geometric phase during the evolution through a 
closed path in the external parameter space |Berry 1984 ISimonj . 



7„(C) =i ^(n,R\V M \n,R) -dR = j A n {R) ■ dR . (301) 

— * 

The path integral in the parameter space, j n (R) , represents the "celebrated" Berry phase 
and the integrand, 

A n (R) =i(n,R\ Vg\n,R) , (302) 

is a vector-valued function called the Berry connection or the Berry vector potential. 

Let us note that the Berry connection transforms as a gauge vector field. Indeed, the 
gauge transformation 

\n,R) — ► exp (i£(R))\n,R) , (303) 
with C,(R) being an arbitrary smooth function, acts on the Berry connection as 

A n (R) — ► A n (R)-V^(R). (304) 

Consequently the Berry vector potential A n (R) transforms as an U(l) gauge potential. 

In analogy to the electrodynamics, the gauge field tensor derived from the Berry vector 
potential reads 

^ = M> A "-^ A '' (305) 

and is known as the Berry curvature. Using Stokes's theorem, we can express the Berry 
phase as an integral of the Berry curvature throughout the surface S enclosed by the path 
C, 

7 " = I ! dR i A dW m . (306) 
2 Js 

It is worth noting that the Berry curvature can be viewed as a U(l) gauge- invariant magnetic 
field in the parameter space. It is therefore observable. 

The non-Abelian gauge potential : degenerate states 

The Berry phase admits a non-Abelian generalization when the energy levels of the Hamil- 
tonian are degenerate [Wilczek 1984j . 

Thus, we now consider a quantum system described by an Hamiltonian TL(R(t)) with 
fc-fold degenerate ground states for all values of the external parameter R(t) . For simplicity, 
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we fix k = 2 such that the energy levels are two-fold degenerate and the Hamiltonian has 
two independent eigenstates, \n a , R(t)) , a = 1, 2 . The eigenvalue equation now reads as 

H(R(t)) K, R(t)) = S n (R(t))\n a , R(t)) , (307) 

where, without loss of generality, the eigenstates are chosen such that 

(n a , R(t)\n b , R(t)) = 5 b , a, b = 1,2. (308) 

Assuming that the system starts in one of its eigenstate \n a , -R(O)) , the adiabatic approx- 
imation implies that the system stays in its initial instantaneous eigenstate, after a cyclic 
tour through the space of parameters. The eigenfunctions of the system are 

\* a m (t)) = \m b ,R(t))UZ(R), (309) 



where U£(R) G SU{2) is an unitary matrix. Let us now substitute the expression (309) 
into the time-dependent Schrodinger equation, 

zh^l(t)) = n(R(t))\^ m (t)) , (310) 

and multiply it from the left by {n a , R(t)\ , one finds 

^A b c (R)U^R)+ i d ^^-^S n U:(R) = 0, 0,6,0=1,2, (311) 

where we introduced the notation, 

A b (R) = i{n c ,R(t)\V n \n b ,R(t)) . (312) 

with a,b,c— 1,2 denoting matrix indices. Hence, the vector potential A(R) is a (2 x 2) 
anti-Hermitian matrix lying in the su(2) Lie algebra. Indeed, under the non-Abelian gauge 
transformation, 

\m' a , R(t)) = \m b , R(t)) U£{R) , (313) 



the field (312) transforms as 



A b (R) — ► U- 1 (l b (R) - U , (314) 

and therefore defines a su(2) valued Berry vector potential. The associated non-Abelian 
Berry curvature then reads as 

F » = m A ^-m^ +i ^ A ^ ■ (315) 
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Writting the solution of (311) in terms of the path-ordered integrals, 

U = Pexp (J A b c {R) ■ dRj x exp J dt' S n {R{t'))\ . (316) 

It is worth remarking that the system undergoes a SU(2) rotation which depends on the 
path taken, 

7n = Pexpf <f l*(R)-dR] . (317) 



The latter defines the non-Abelian generalization of the Berry phase factor known as the 
Wilson loop. 

Compared to the Berry phase which is always associated with a closed path, the Berry 
curvature is truly a local quantity. It provides a local description of the geometric properties 
of the parameter space while the Berry phase can be identified with the holonomy of the 
fiber bundle [Si mon] . 

Also, the Berry curvature also plays the role of a (non-Abelian) magnetic-like field, 
which affects the particle dynamics in his neighborhood. A relevant example is provided 
with the "Berry" non-Abelian monopole-like fields arising in diatomic molecule systems 
[Wilczek 1986] . see the next section. 
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C. Monopole-like fields in the diatomic molecule 



As first investigated by Moody, Shapere and Wilczek [Wilczek 1986] . a "truly" non- 
Abelian gauge fields mimicking monopole-like fields can arise in a diatomic molecule system. 
These come from the non-Abelian generalization |WilczekT9 84] of the Berry gauge poten- 
tials. In this case, we consider sets of levels, /c-fold degenerate, subjected to adiabatically 
varying external parameters. For k — 1, a single level, we recover the U(l) gauge fields 
discussed by Berry and Simon |Berry 1984 ISimon] . For k > 2 , the effective gauge fields 
take a "truly" non-Abelian form. 

The latter can be extended to systems where the slow dynamical variables are no longer 
external but are themselves quantized. This is the case, in particular, for the diatomic 
molecule where the quantized parameters define nuclear coordinates [Wilcze k 1986J. 

To see this, let us consider a diatom with two atomic nuclei and one or more gravitating 
electrons. The study of this system amounts to investigating a many-body problem which 
reduces, in the simplest case, to a three-body problem. Neglecting the spin degree of freedom 
and the relativistic effects, the Hamiltonian employed in the diatomic molecule system reads 
[in units h — 1 ] , 

f 1 n 2 

H(Xi, x k ) 

I LIU l. ' ' L . \ I , I 

(318) 

3=1 k>j 

Here the atomic number Zj corresponds to the electric charge of the jth nucleus and the 
positions Xj and x& denote the nuclei and the electrons coordinates, respectively. The 
bracketed terms in the Hamiltonian are the kinetic energy of the electron of mass m plus 
the kinetic energy of the nuclei of mass Mj , with V Xi and Vx fe referring to the Laplacians 




of the ith electron and of the kth nucleus, respectively. The two following terms in (318) 
define the classical Coulomb electron-nuclei interaction and the nucleus-nucleus interaction, 
respectively. The remaining term represents the electron-electron interactions. 

From now on we consider the simple configuration of the molecular ion which pos- 
sesses only one gravitating electron and two identical hydrogen nuclei. Then, the electric 
charge of the nuclei are the same, Z\ — Z2 — Z , and electron-electron interactions van- 
ish since only one electron is considered in the present context. The total non-relativistic 
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wave function ^(x, Xi) of this diatomic molecule system is a solution of the stationary 
Schrodinger equation, 

U Xi) = S Xi) . (319) 

The description of the diatomic molecule properties is commonly made using the Born- 
Oppenheimer approximation. Indeed, the Born-Oppenheimer or adiabatic approximation 
is applied to separate, in an appropriate way, the electronic motion and the slower nuclei's 
degrees of freedom that couple to it, since Mj ^> m . To investigate electronic motions, we 
first assume that the nuclei positions X\ and X2 are fixed and correspond to infinite nuclear 
masses Mi = M 2 = 00 . Thus for a fixed nuclear configuration, we obtain the electronic 
Hamiltonian, H e i , carrying a parametrical dependence on the nuclear relative coordinate 
X 12 — Xi — X 2 , 



n el (x,x 12 ) = _-Lvi + v(£,Xi 2 ) , 

V(x,X 12 ) = -— - — + 



(320) 



\X 1 -x\ \X 2 -x\ \X 12 \ 
Since X\ 2 is just a parameter, then the last term in the previous potential is a constant 
and shifts the eigenvalues only by some constant amount. In the context of the Born- 
Oppenheimer approximation, we consider also 

*(£, Xi) » X 12 ) , (321) 

where the molecular wave function, fy(x, X i2 ) , can be expanded into a combination of the 
electronic wave function ip m (x) and the nuclear wave function Xm(Xi 2 ) , 

X 12 ) = <Pm& X 12 ) X m(X 12 ) ■ (322) 

m 

Note that the electronic eigenfunction depends implicitly on the nuclear relative coordinate, 
X 12 , and the summation index m denote the eventual energy's degeneracy of the electronic 
eigenstate. Hence, the electronic eigenfunction obeys to the electronic stationary Schrodinger 
equation, 

"He/ <fm(x) = £el,m ^>rn{x) , (323) 

and form a complete set. While when investigating the nuclear motions, we have to consider 
the electron as remaining in the same quantum eigenstate so that the nuclear wave function 
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is a solution of the Schrodinger equation with an effective potential generated by the electron, 

E X k(X 12 ). (324) 



Let us first investigate the eigenvalues equation (323) for the molecular ion H 2 with the 



nuclei located on the orthogonal z -axis [see Figure [7] . We use the spherical coordinates 



(r, 9, (f) [see Figure [Tj to rewrite the electronic Hamiltonian (320) in this coordinate system, 
n a (r,9,. 

V(r,6,(/)) = v(r,e) = 



r H V V (r, 9, 

2m r dr 2 2m r 2 

Ze 



Ze 



^r 2 + R 2 2 + 2rR 2 cos 9 ^r 2 + R 2 -2rR l cos 9 




Electron 



■■■<!> 



FIG. 7: Molecular ion H 2 with the set (r, 9, <fi) representing spherical coordinates and 
{9 a , 6b, Ri, R2, r a , rt,} providing us with elliptic coordinate system on the plan. 



The Casimir L 2 and the projection of the electronic orbital angular momentum L z read 



Lr 



d 2 



1 d 



1 d 2 



d 



(325) 



d9 2 tan 9 89 sin 2 9 dcf) 2 ' z d<p ' 

Note that the potential, V(r, 9) , which does not depend on the azimuthal angle , is 
rotationally symmetric around the axis of the nuclei, 



[L z , V(r, 







(326) 
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Since the component of the angular momentum L z also satisfies, 

[L z , L 2 ] = 0, 



then, (326) and (327) yields the conservation of L z along the electronic motion, 

[L z ,H a ] = 0. 



(327) 



(328) 



The eigenvalues of the quantized quantity L z are given by the eigenvalues equation, 

L z ip m (r, 9, 0) = mip m (r, 9, 0) . (329) 

The diatomic molecule therefore possesses a privileged direction carried by the axis of the 
nuclei [the z-axis]. Thus, L z , generates an 5*0(2) symmetry group. Moreover the config- 
uration of the nuclei is also invariant under spatial inversion. Consequently the electronic 
Hamiltonian, Tiei , admits the same symmetry group Q as the nuclei's configuration, 



G = 50(2) x (Parity) . 



(330) 



Taking into account the 5*0(2) symmetry of the diatomic molecule, we can separate the 
electronic wave functions, (p m , under the form of the product, 



<Pm{r, 0, 0) = g m ((f>)f(r, 9) . 



(331) 



Injecting the expanded form (331) into the eigenvalues equation (329), it is straightforward 



to obtain the normalized eigenfunction g m , 

1 



9±m (0) 



exp (±z m(j 



m e Z . 



(332) 



Now, our task is to investigate the function, f(r, 9^j, appearing as a part of the electronic 



wave function <p m (r, 9, 0) in (331). To this, we switch to elliptic coordinates, (r, i 
(f, 7]), see Figure |7| 

p _ r_g + n _ r a - r h 
K ~ R ' V ~ R ' 

R = R 2 + R U £ G [1, oo[ , 77 E [-1, 1] . 

We can now express the potential and the Laplacian operator in elliptic coordinates, 

4eZ £ 



(333) 



R £ 2 -?7 2 



R 2 (£ 2 - if 



{(e-i) 



d 2 d 2 

Q£2 K I > drj 2 



(334) 



1 



1 



d 2 



+2£--2r^4 , , 

<9£ dr] \£ 2 — 1 7] 2 — 1 / d(j) 2 
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so that the electronic Schrodinger equation (323) which takes, in elliptic coordinates, the 
form 

Uei g± m (0) / (f , rj) = E e i 9± m (0) / (£, v) , (335) 

separates into 



mR? 
2 



((" 2 - i »l? +2 " 



<9 



m 



(336) 



mi? 2 V <9?7 2 drj 77 2 — 1 

Let us note that the dependence on the azimuthal angle has disappeared from the bracketed 



terms in (336) and is replaced by the parameter m 2 . Consequently, the electronic energy 
spectrum depend on m 2 , 



el,±m ^ei,=Fm j 



(337) 



so that, for m/0, each electronic level is doubly degenerated. Moreover, it is now clear 
that the variables £ and rj separate in the eigenfunctions /(£, rf) as 



/ (£, r?) = /o (0 A (V) 



(338) 



where /o and f\ are solutions of the following spheroidal wave equations 6 , 



(e - 1) 

(v 2 - 1) 



d 2 



d 



-2 t-2 



m 



d 



drf ' ^ c??7 



2 2 

-a + p r\ — 



m 



rf — 1 



/o(0 = o 

A(»?) = o. 



Here m 2 and a are the separation constants of the differential equation, with 

R 2 

^ = 2Rme 2 , p 2 = -—mS e i. 



(339) 



(340) 



The bound states are solutions of the Schrodinger equation (323) associated with quantized 
negative energies, {E eliTn } meM . 



Collecting our results (332) and (339) provides us with the complete electronic wave 
functions, 



1 



fo (0 h (v) ex P (±im<i>) , ra G Z . 



(341) 



Spheroidal wave equations are generalization of Mathieu differential equations. 
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Let us recall that the integer m which corresponds to the eigenvalue of the component 
of the angular momentum along the axis of symmetry z is a "good" quantum number, 



L z ip m = ±A ip m , with A = \m\. (342) 

When introducing the electronic spin degree of freedom, S z , a non-vanishing additional 
term can be included, at each nuclear configuration, in the electronic Hamiltonian of the 
diatom , namely 

H so = /iL z ■ S z , fjt e R, (343) 

corresponding to the spin-orbit effects. In that event, we can show that S z and L z are 
separately conserved, 

[5*2, H e i\ = [S z , L z ■ S z ] — 

(344) 

[L z , Uei] = [L z , L z -S z ] = 0. 
The projection of the total angular momentum, J z = L z + S z , onto the z -axis is therefore 
quantized, 

Jz<P±k(€, V, <t>) = k<p±k(£,, V, <t>) , k = ±A, (345) 

so that for fixed quantum number Ao , the eigenvalue, A , associated with J z takes the two 
half-integer values, 

A = Ao-^, A + ^. (346) 

Consequently, the introduction of the electron spin degree of freedom does not modify the 
double degeneracy of the electron system. The novelty here is that even the ground level 
Ao = remains doubly degenerated. Thus, the electronic wave functions are now charac- 
terized by the quantum number A , 

V±k (f , V, <P) = ((£, V, 0) I ± k) = -4=/o (0 h (V) exp (±t k<j>), k = ±A . (347) 

V ^7T 



Following the Born-Oppenheimer approximation, after describing the electronic wave func- 
tions labeled by the index k with the energy eigenvalues S e ^k which are parametric func- 
tion of the relative internu clear coordinate, we are now interested on the nuclear motions 



described by the Schrodinger equation (324). Let us sandwich (324) between electronic 
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eigenstates, 

f 2 1 Z 2 

~ / dx(p*(x, X 12 ) ( ^ WTJ^X, ~ Kel - -f—) ^ X X2 )Xm{X 12 ) 

J i= i ZM \Xu\ m (348) 

J m 

Thus, this provides us with 

2 1 f 

~^2ohf dx^ n (x,X 12 )V 2 Xi Y^ m (x,X 12 )x m (X 12 ) 

^ J z2 ™ (349) 

+£ dX n(X 12 ) + -^—Xn{X 12 ) = E X n{X l2 ) , 

\Xi 2 \ 

where the first term can be expanded using the Leibniz rule on differentiation so that we 
obtain for \<p m ) = \m) , 

1 Z 2 
~ (Z)< n l V L' m ^ m + 2 J2< n \ V xJrn)Vxm + % 12 Xn) + fa + 7^~)Xn = E X n ■ 

Hence, we derive the effective Hamiltonian describing the nuclear motions, 

^™ = -^E fafc, + < n I ^x 12 1 k » (V^ 2 + ( k | V^ i2 |m» + (-|- + E d )5 nm . 

In the adiabatic approximation, where the nuclei move slowly when compare to the electronic 
motion, the electron has to be considered to remain in the same 2 -fold degenerate nth level. 
Consequently, the off-diagonal transition terms are neglected, implying the relevant effective 
nuclear Hamiltonian, 

z 2 ( 35 °) 

with A = i( n | Vj£ l2 | n ) and V(Xi 2 ) = — V £ e i ■ 

\X\ 2 \ 

Here V acts as an effective scalar potential for nuclear motion and the induced gauge 
potential A is a (2 x 2) matrix, since the state | n ) belongs to a 2 -fold degenerate level, 



see (337) and (346). Hence, A transforms as a £7(2) gauge potential. 

For the nuclear axis in the initial direction given by the polar and azimuthal angles 
= = 0, 

\n(e z )) = |n(0,0)>, (351) 
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we can generate a set of eigenstates adapted to nuclei pointing toward (9, 0) by rotating the 
initial eigenstate. Then, the Wigner theorem provides us with the two possible parametriza- 
tions, 

\n(9, 0)) = exp(iJ 3 0) exp(iJ 1 6 l ) exp(— iJ 3 (fr)\n(0, 0)), for 9 ^ n 



(352) 



\n(9, 0)) = exp(iJ 3 0) exp(iJi^) exp(iJ 3 0)|n(O, 0)), for 9 ^ . 
Note that the two previous parametrizations are linked by 

| n) — exp(2in0)| n) . 



(353) 



Thus, the U (2) gauge potentials, which are defined on the space spanned by the electronic 
eigenstate, depend on the geometry of the 2 -fold degenerate eigenstate space so that 

A r = i{ n(r, <j>)\d r \ n(r, 0, </>)) 

A<, = i{n(r,e,<i>)\de\n(r,e,<l>)) ( 354 ) 
= i{n(r,9,(f>) \ <9 \ n(r, 0,0)) . 
Performing the calculation in the case of the 9 ^ n parametrization 7 leads to 
A r = 0, 

A g = (n(0,0) | -cos0Ji + sin0J 2 |n(O,O)), (355) 
A 4> = (n(0,0) | (1 -cos 9)J 3 + sin (sin 0Ji + cos 0J 2 ) | n (0, 0) ) . 



We posit, 



Ji = aa u J 2 = ba 2 , J 3 = ca 3 , 



(356) 



and we obtain by direct calculation the shape of the non-vanishing gauge potentials induced 
by nuclear rotations, 



Ag = —a cos <j\ + b sin cr 2 , 

A^ = sin 9 (a sin + 6 cos <x 2 ) + c (1 — cos 9) a 3 . 
The corresponding field strength, Fg^ , reads 

Fq^ = a sin # a 3 + (cos 6* — 1 ) (/5 cos <r 2 + 7 sin o"i) , 



(357) 



(358) 



7 The procedure is exactly the same for the parametrization with 8^0. 



81 



where a, (3 and 7 satisfy the relations, 



a = c — 2ab , (3 = b — 2ac , 7 = — 26c . 



(359) 



Let us now inquire about the real nature of the U(2) gauge potentials (357) induced by 
nuclear motions. Are these imbedded Abelian gauge fields into U (2) ? or not? To respond 
to this question, let us recall that a field strength, 



P 



ma x + na 2 + pa 3 



(360) 



can always be gauge-transformed so that it points in one single direction <j\ , a 2 or cr 3 say. 
In the present context, we search for gauge transformations which rotate the field strength 



(360) in the "Abelian" direction 03 



Pe,4> = jfaa , / ^ 



(361) 



In the limiting case where / = , i.e. in the null field strength configuration, the gauge 
potentials are pure gauge. Then, the gauge potentials must be gauge equivalent to that of 
the vacuum. 

We are looking for matrices U taking values in £7(2), 



U 




(362) 



so that J 7 ^ = UF'q^ U~ 
be solved 



Consequently, we derive from (IVC) the series of constraints to 



(p-f)A + 


(m + in) B 


= 


(m — in) A - 


-(f+p)B 


= 


(m — in) C - 


Hf~P)D 


= 


k (f+P)C + 


(m + in) D 


= 



(363) 



The constraints (S) can be solved provided its determinant vanished, 

det(S) = <^=> f = m 2 + n 2 + p 2 . 



(364) 



We then obtain an equivalence between the length of Tq$ and J 7 ^ which express as con- 
servation of the length of the field strength under a gauge transformation . Solving the 
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constraints (S) , for / ^ , yields 



17 



7 + p 

2/ 



m - in , v \ 

exp (ut) — exp u/il 

/ + p 

V~7^]T exp H ex P (2z/) y 



with arbitrary real constants 



/x = argfA) , i/ = arg(D) 



Applying the previous gauge transformation to (358) and (357) we must take 



m = 7 sin (cos 6* — l) 
n = / 9cos^(cos6 l — l) 
p = a sin 9 , 

v. 

so that the length of the field strength reads 



a 



2 sin 2 9 + ( cos 6 - 1) ( 7 2 sin 2 + /3 2 cos 2 



Without loss of generality, we can choose [Wilczek 1986| IRho 19 92] , 

a = ^(l-K 2 ), /3 = 7 = 8 with kgR 9 , 



so that applying the gauge transformation on the Berry potentials (357) 

A e = U{A e + id e )U- 1 

A <p = U(A (j> + id (j) )U- 1 , 
provides us with the gauge-equivalent potentials, 



A e = T 1 







and 



A„ 



I -(l-cos0) ^-\K\sm6 e^-^^ 



±-\n\smde i{u -^- <t>) 



1 



V^2 



1 — COS I 



(365) 



(366) 



(367) 



(368) 



(369) 



(370) 



(371) 



8 This choice implies that a = b — (+|Ji| — ) = ±^ and c = h . 

9 Since the electronic eigenstates are not eigenfunctions of angular momentum, but only of J3 , k can take 
any real value. 
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It is now clear that the Berry gauge potentials (357) or (370) and (371 ) become "Abelianized" 
gauge potential for k = . In that event, they represent a Dirac monopole field of unit charge 
imbedded into U{2) , 



A ft = 



A d 



1 



1 — cos6) 03 . 



(372) 



For \k\ , we obtain the truly non-Abelian case, where the off-diagonal terms can not be 
eliminated 10 , 



A e = =p— ( cos <p a-i - sin <\> a 2 ) , 



1 



(373) 



A^ = ±— sin^(sin0cri + cos0o"2) + r (l — cos 0)03 . 



The corresponding field strength is 



(l — k 2 ) sin #(73 . 



(374) 



The field strength (374) superficially resembles to that of a monopole field but the interpre- 
tation is quite different. Indeed, \k\ ^ is not quantized here and the gauge fields induced 
by nuclear motions of the diatomic molecule are truly non-Abelian [Wilczek 1986] . See also 
|Zygelman 1990 . 



Note that when k — ±1, the field strength vanishes and (373) is a gauge transform of 
the vacuum. 



Our next step is to present the monopole-like field (373) in a more convenient "hedgehog" 
form. This can be achieved, by applying a suitable gauge transformations |Jackiw 1986J to 



the diatomic molecule gauge potential (373). Finally, the Berry gauge potential mimics the 



structure of a non-Abelian monopole |'t Hooft 19741 Polyakov 1974 



v l K)tiaj ^ , 



(1 - K 2 )e ijk — 



(375) 



Note the presence of the unquantized constant factor (1 — k 2 ) in the above magnetic field. 



Classical dynamics and conserved quantities 



Now we turn to investigating the symmetries of an isospin-carrying particle, with unit 



charge, evolving in the monopole-like field of the diatom (375) plus a scalar potential. The 



10 



Here we fixed fj, = v . 
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Hamiltonian describing the dynamics of this particle is expressed as 

n = Y~ ^ /4 ) €ijk ^ ^ + K(f ' X<1) ' ^=Vi~m\ (376) 

where the spin-rotation coupling disappears when we study particle carrying null gyromag- 
netic ratio, g = . The resulting Hamiltonian has the same form of that of a scalar particle 
11 evolving in the same magnetic field. We define the covariant Poisson-brackets as 

dN dM ~ DM dN DM dN 

{M,N} = Dj M Wj - _ DjW + 2.F J J__ (377) 

where Dj is the covariant derivative, 

Drf = dif - e abc X a A) ?L . (378) 
The commutator of the covariant derivatives is recorded as 

[A, D s ] = -e abc l a F^ . (379) 

The non-vanishing brackets are 

{x\7r J } = 5}, {7r J ,7r J }=X a ^., {l a ,l b } = -e abc l c , (380) 
and the equations of motion governing an isospin-carrying particle in the static non-Abelian 



gauge field ( |375[ ) read 

' Xi - X a F^ x j + DiV = 
± a + e afec Z b fi^-|^) =0 



(381) 



The first equation in (381 ) describes the 3D real motion implying a generalized Lorentz force 
plus an interaction with the scalar potential; while the second equation is the Kerner-Wong 
equation augmented with a scalar field interaction. The latter describes, as expected, the 
isospin classical motion. 

Let us now recall the van Holten procedure yielding the conserved quantities. The con- 
stants of the motion are expanded in powers of the momenta, 

Q(f,X,vf) = C(x,l) + Ciix,!)^ + -CijixfiniKj + ■■■ , (382) 



i.e. particle without spin. 
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and we require Q to Poisson-commute with the Hamiltonian, 



n 2 



{Q, n = Y + v(x,r)} = o. (383) 

We therefore get the set of constraints which have to be solved, 

BC BV 

W + e^-^O, o(0) 

.da dV 



D,C = l a F°jCj + C V D S V + e abc l a — \ — , o(l) 



~ ~ dC dV 

DiCj + Djd = X a {F1 k C kj + Fj k C k i) + C ijk D k V + 6 <>^ a ~g^"Q^> °( 2 ) (334) 



DiCjk + DjCki + DtCij — Z a (F%Cijk + F^Ciki + F&Cuj 

dCi jk d\ 



+C ijM D l V + e abc l" C ^^, o(3) 



Turning to the zeroth-order conserved charge, we note that, for k 7^ 0, the used-to-be 
electric charge, 

— * 

— » *r 

Q = —, (385) 
r 

is not more covariantly conserved in general, 

{Q,H}=n- DQ, DjQ = - (l j - Q^j . (386) 



An exception occurs when the isospin is aligned into the radial direction, as seen from (386). 
A detailed calculation shows that the equation DjQ = can only be solved, for imbedded 
Abelian monopole field, when k — 0, ±1 . 
Nor is Q 2 conserved, 

{Q 2 ,n} =2kQ(tt-DQ). (387) 
Note for further reference that, unlike Q 2 , the length of the isospin, X 2 , is conserved, 

{H,X 2 } = 0. 



The monopole-like gauge field (375) is rotationally symmetric and an isosp in-carrying 
particle moving in it admits a conserved angular momentum |Wilczek 1986| IJackiw 19 86] . 
Its form is, however, somewhat unconventional, and we re-derive it, therefore, in detail 
Nsrome 02/20091 . 
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1) We start our investigation with conserved quantities which are linear in the covariant 
momentum. We have therefore 



(388) 



so that the series of constraints (384) reduce to 

QD i V + e abc l a 



dX b dl c 



0, 



o(0) 



~ 8(1 f)V 

A^I^. + e^-, o(l) 

k DiCj + Djd = . o(2) 
We use the Killing vector generating spatial rotations, 

C = n x x . 



(389) 



(390) 



Choosing V = V(r), we see that, again due to the non-conservation of Q , DjV ^ in 



general. The zeroth-order condition C ■ DV = in (389) is, nevertheless, satisfied when 



V is a radial function independent of X , since then DV = W , which is perpendicular to 
infinitesimal rotations, C. 



Evaluating the right hand side of the first-order constraint of (389) with F^ k as given in 



(375), the equation to be solved becomes 

DiC = (1 - K 2 )- ( (n--)— -n, 



(391) 



a; 



In the Wu-Yang case, k = 0, this equation was solved by C = —n ■ Q— . But for k ^ , 

r 



the electric charge, Q , is not conserved, and using (386), (391), as well as the relations 

Q 



Di(X-n) 



A Qn 



x 



(1 — K 

Q 



n ■ r 



Hi 



rii - (1 + /t)(n • r)— J + -(n • f)Xj 

;i-« 2 )^, 



(392) 



we find, 



;i-«)a gn 



KDiiX-n) +DiC. 



This allows us to infer that 



C7= -((l-«)Q^ + /«z) -n. 



(393) 
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The conserved angular momentum is, therefore, 



X X 7T 



'1 - «) Q - + «T =Q- + k 
r r 



x -A x 
-xx x — 
r I r 



(394) 



consistently with the results in |Jackiw 1986} IRho 1992] . Moody, Shapere and Wilczek 
[Wilczek 1986] found the correct expression, (394), for k = but, as they say it, "they 
are not aware of a canonical derivation when k ^ 0". Our construction here is an alter- 
native to that of Jackiw [Jackiw 1986j . who obtained it using the method of Reference 
[Jackiw-Manton 1980] . In his approach, based on the study of symmetric gauge fields 



Forgacs-Manton 1980 , each infinitesimal rotation, (390), is a symmetry of the monopole in 



the sense that it changes the potential by a surface term. 

It is worth noting that comparison with the Wu-Yang case yields the "replacement rule" , 



x 



r 



(395) 



For k = we recover the Wu-Yang expression (283). Eliminating tt in favor of p — A = tx 



allows us to rewrite the total angular momentum as 



J = x x p — Z. 



making manifest the celebrated "spin from isospin term" [Jackiw 1976] . 



(396) 



Alternatively, a direct calculation, using the same formulae (386)-(392), allows us to 
confirm that J commutes with the Hamiltonian, {J^'H} = 0. 



x 



Multiplying (394) by — yields, once again, the relation (290) i.e., 

r 



J ■ 



x 



-Q 



(397) 



the same as in the Wu-Yang case. This is, however, less useful as before, since Q is not a 
constant of the motion so that the angle between J and the radius vector, x(t), is not more 



a constant. The components of the angular momentum (394) close, nevertheless, to so(3) 



In addition of J , it is worth mentioning that the Casimir 

J 2 = (x x vf) 2 + (1 - k) 2 Q 2 - K 2 ! 2 - 2k J- X 
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(398) 



(399) 



is obviously conserved since the angular momentum of the diatom is conserved ( 394 ) . 



2) Returning to the van Holten algorithm, quadratic conserved quantities are sought by 
taking 

C ijk = C m ■ ■ ■ = . (400) 



Consequently (384) reduces to 



DiC = T a F«f 3 + CuDiV + e abc I 



o(0) 
o(l) 



~ ~ BC- ■ BV 

AC, + Did = l a (F? k C kj + Ff k C ki ) + e abc X a ^ — , o(2) 



(401) 



DiCjk + DjCki + D^Cij — . 
We consider the rank- 2 Killing tensor, 



o(3) 



C^j %2 — ^ &3 ^ ^ 2 



(402) 



which satisfies the third-order constraint of (401). Injecting (402) into the second-order 
constraint yields, 

DiCj + DjQ^O, (403) 
which can be solved by taking Cj = . For radial potentials independent of X, it is straight- 



forward to satisfy the first- and the zeroth-order constraints of (401) with C — 0. Thus, we 
obtain the conserved Casimir, 



(x x 7?) 2 = x 2 ti 2 — (x ■ 7?)' 



(404) 



which is the square of the non-conserved orbital angular momentum, L = x x tt . 

Since J 2 and L 2 are both conserved, it is now straightforward to identify the charge, 



r = J 2 - L 2 = (l - k) 2 Q 2 - k 2 ! 2 - 2kJ- X 



(405) 



which is conserved along the motion in the monopole-like field of diatomic molecule. It is 
worth noting that the charge V corresponds, in the Abelian limit with k = , to the square 
of the electric charge. As the constants of the motion J , J 2 and L 2 , the charge T is 
conserved for any radially symmetric potential, V(r) . 
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Note that T can also be obtained by using the Killing vector, 

C = 2k(xxX), (406) 



into the van Holten algorithm (389). 



Let us now decompose the covariant momentum, into radial and transverse components, 
with the vector identity, 

L 2 



(tT) 2 = (7f • If + (7? X = 7T r 2 + ^ . (407) 



This hence allows us to express the diatomic molecule Hamiltonian (376) as 



n = i<* ■ f )' + £ - { ('-^-^- 2kJ .x | + nr) (408) 

Suggesting that the charge takes the fixed value Q 2 = X 2 = 1/4, Jackiw found a similar 



decomposition as (408) [Jackiw 19 86J. but this is, however, only legitimate when k — 0, 
since Q 2 is not conserved for k^O. 

For k^O, the "good" approach is to recognize the fixed charge T , which yields the nice 
decomposition, 

« = 5(*-f) > + ^-^ + W- < 409 > 

Let us underline that the effective field of a diatomic molecule provides us with an in- 
teresting generalization of the Wu-Yang monopole. For k ^ 0, ±1, it is truly non-Abelian, 
i.e., not reducible to one on an U(l) bundle. No covariantly constant direction field, and, 
therefore, no conserved electric charge does exist in this case. 



The field is nevertheless radially symmetric, but the conserved angular momentum (394) 
has a non-conventional form. 

In bundle terms, the action of a symmetry generator can be lifted to the bundle so that it 
preserves the connection form which represents the potential. But the group structure may 
not be conserved; this requires another, consistency condition [Jackiw-Manton 1980] , which 
may or may not be satisfied. In the diatomic case, it is not satisfied when k^0, ±1. 

Is it possible to redefine the "lift" so that the group structure be preserved ? In the 
Abelian case, the answer can be given in cohomological terms [Duval 1982j . If this obstruc- 
tion does not vanish, it is only a central extension that acts on the bundle. 

In the truly non-Abelian case, the consistency condition involves the covariant, rather 
than ordinary derivative and covariantly constant sections only exist in exceptional cases - 
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namely when the bundle is reducible. Thus, only some (non-central extension) acts on the 
bundle. 



It is worth noting that for k 7^ the configuration (375 ) does not satisfy the vacuum Yang- 



Mills equations. It only satisfies indeed with a suitable conserved current |Jackiw 1 986]. 

. k(1-k 2 ^ 



xxT 



i-^ik Jki J 



Interestingly, this current can also be produced by a hedgehog Higgs field, 



(410) 



Jk 



vT 



K X. 



a 



(411) 



For k — , it is straightforward to derive the conserved Runge-Lenz vector since this case 
is exactly equivalent to the Wu-Yang case, an imbedded Abelian monopole. For k 7^ , ±1, 



we derived a new conserved charge, namely V , which has an unconventional form, see (405). 
In the limit case k = , this conserved charge reduces to T = Q 2 ; while for k = ±1 , we 
obtain V ~ L ■ X . 



Let us emphasize that the derivation of the non- Abelian field configuration (375) from 
molecular physics [W ilczek 1986] indicates that our analysis may not be of purely aca- 
demic interest. The situation could well be analogous to what happened before with the 
non- Abelian Aharonov-Bohm experiment, first put forward and studied theoretically in 
|Wu Yang 1975 , Horvathy 04/198"5] , but which became recently accessible experimentally, 
namely by applying laser beams to cold atoms Ohberg 2005 Ohberg 2007 IDalibard et "aT) . 
A similar technique can be used to create non- Abelian monopole-type fields [Dalibard et al.\ . 



V. SUPERSYMMETRIC EXTENSION OF THE VAN HOLTEN ALGORITHM 

We investigate the super and dynamical symmetries of a fermion in external magnetic 
fields using a SUSY extension of the van Holten framework, based on Grassmann-valued 
Killing tensors. 



A. Supersymmetry of the monopole 

In this section, we investigate the super- and the dynamical symmetries of fermions in a 
D -dimensional monopole background. Following an interesting result of D'Hoker and Vinet 
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[D'Hoker 1984] . a non-relativistic spin-| charged particle with gyromagnetic ratio g = 2 
interacting with a point magnetic monopole, admits an osp(l|2) supersymmetry. This was 
also seen in the following papers [Gibbons 19931 DeJonghe 1995 Rietdijk Horvathy 2000 



Plyushchay 20001 |Plyushchay 04/2000} ILeiva 2003 1 |Horvathy 20051 - 



Later, Feher |Feher 1987] has shown that a g = 2 spin-particle in a monopole field does 
not admit a Runge-Lenz type dynamical symmetry. 



Another, surprising, result of D'Hoker and Vinet D'Hoker 01/1985, D'Hoker 09/1985 



D'Hoker 04/T986] says, however, that a non-relativistic spin-| charged particle with anoma- 
lous gyromagnetic ratio g = 4 , interacting with a point magnetic monopole plus a Coulomb 
plus a fine-tuned inverse-square potential, does have such a dynamical symmetry. This is to 
be compared with the one about the 0(4) symmetry of a scalar particle in such a combined 
field [Mcintosh 19701 Zwanziger 1968| . Replacing the scalar particle by a spin 1/2 particle 
with gyromagnetic ratio g = 0, one can prove that two anomalous systems, the one with 
g — 4 and the one with g = are, in fact, superpartners [Feher 1988] . Note that for both 
particular g-values, one also has an additional o(3) "spin" symmetry. 

On the other hand, it has been shown by Spector |Spector| that the M = 1 supersym- 
metry only allows g — 2 and no scalar potential. Runge-Lenz and SUSY appear, hence, 
inconsistent. 

We study the bosonic as well as supersymmetries of the Pauli-type Hamiltonian, 



5! 



eg 



S-B + V(r), U = p-eA 



(412) 



which describes the motion of a fermion with spin S and electric charge e , in the combined 
magnetic field, B, plus a spherically symmetric scalar field V(r), which also includes a 



Coulomb term (a "dyon" in what follows). In the Hamiltonian (412), II denotes the gauge 
covariant momentum and the constant parameter g represents the gyromagnetic ratio of 
the spinning particle. 



Let us first describe the Hamiltonian dynamics, defined by (412), of the charged spin-^ 
particle, moving in the flat manifold Ai D+d . Note that Ai D+d is the extension of the 
bosonic configuration space Ai D by a d -dimensional internal space carrying the fermionic 
degrees of freedom |Cariglial . The (D + <i)-dimensional space Ai D+d is described by the 
local coordinates (x M , ip a ) where fi = ,D and a = 1, • • • ,d. The motion of the 



92 



spin-particle is, therefore, described by the curve 

t->(x(t), V(r)) G M D+d . (413) 

We choose D = d = 3 and we focus our attention to the spin-| charged particle inter- 
acting with the static U(l) monopole background, 

B = VxA = ^, (414) 



so that the system is defined by the Hamiltonian (412). We introduce the covariant hamilto- 
nian formalism extending van Holten's framework to fermions. The basic phase-space reads 
(x J , Uj, ip a ) , where the variables tp a transform as tangent vectors and satisfy the Grassmann 
algebra, 

^ip j + V> j y = . (415) 

The internal angular momentum of the particle can also be described in terms of vector-like 
Grassmann variables, 



S* = -£«W- (416) 
Defining the covariant Poisson-Dirac brackets for functions / and h of the phase-space as 

wS-S^^TO^- 1 ^' (417) 

where a* = (0, 1) is the Grassmann parity of the phase-space function / and the magnetic 
field reads = (l/2)eyn,i^-fc . It is straightforward to obtain the non- vanishing fundamental 
brackets, 

{x\ ILj = 5), {IL, n,} = eF iv {^, ft} = , (418) 

[S\ G 1 } = e k ij G k with G k = *fj k , S k . (419) 

It follows that, away from the monopole's location, the Jacobi identities are verified 
|Jackiw 1985t IChaichian 20 09J. Thus, the equations of motion can be obtained in this co- 
variant Hamiltonian framework 12 , 

G=|Gx5, (420) 
n = ellx B- W(r) + y v(S-B} . (421) 



12 



The dot means derivative w.r.t. the evolution parameter, 4- . 
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Equation (420) shows that the fermionic vectors S and ip are conserved when the spin and 



the magnetic field are uncoupled, i.e. for vanishing gyromagnetic ratio, g = . Note that, 



in addition to the magnetic field term, the Lorentz equation (421) also involves a potential 



term augmented with a spin-field interaction term (Stern and Gerlach term). 

We now proceed by deducing, in a classical framework, the supersymmetries and con- 



servation laws of the system (412), using the SUSY extension of the van Holten algorithm 



Ngome 03/2010| developed in section |IIB| What is new here is that the generators of SUSY 



are Grassmann-valued Killing tensors. We expand the phase-space function, associated with 
one (super) symmetry, in powers of the covariant momenta, 



n— 1 

q (x, n, = c(x, $ + X) m ctl "' tk & $ n ^ • • • n ^ • 

fe=l 



(422) 



Note the dependence on Grassmann variables of the tensors C(x, if)) . Requiring that Q 
Poisson-commutes with the Hamiltonian, {"H s , Q] = , implies the series of constraints, 



CAV + l ^ l i) m c 3 d 3 F lm - e ^ m 



0. 



djC = C jk d k V + eF jk C k + l -^L^ m C jk d k F lm 



dijj a 



^(jCfc) C j km d m V -\- c (Fj m C mk -\- F km C l 
+ % ^i^ c .. kd . Fir 



mjj 



eg dCj, 



2 T d^ a 

9(j Cfci) Cj k i m O m V -\- c {Fj m C mk i -\- Fi m C m j k -\- F km C m ij^ 



o(0) 
o(l) 

o(2) 
o(3) 



(423) 



This series of constraint can be truncated at a finite order n provided the higher order 
constraint becomes a Killing equation. The zeroth-order equation can be interpreted as 
a consistency condition between the potential and the (super)invariant. Apart from the 
zeroth-order constants of the motion, i.e., such that do not depend on the momentum, all 



other order-n (super)invariants are deduced by the systematic method (423) implying rank 



n Killing tensors. Each Killing tensor solves the higher order constraint of (423) and can 
therefore generate a (super) invariant. 

We focus our attention on searching for conserved quantities which are linear or quadratic 
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in the covariant momenta. Thus, we have to determine generic Grassmann-valued Killing 
tensors of rank-one and rank-two. 

• Let us first consider the Killing equation, 



(424) 



Following Berezin and Marinov [Berezinj . any tensor which takes its values in the Grassmann 
algebra may be represented as a finite sum of homogeneous monomials, 



(425) 



k>0 



where the coefficients tensors, C\ a , are completely anti-symmetric in the fermionic indices 



{afc} . The tensors (425) satisfy (424), from which we deduce that their (tensor) coefficients 
satisfy 



d t d,C k ai ... a Jx) = 



(d j e h av .. am {x) + d k ci v .. am (x)) r i ---r- = o = 

providing us with the most general rank-1 Grassmann-valued Killing tensor 
C\x, $) = { Mij x j + N*) ai ... a ip ai ■■■^ a \ M ij = -M ji 



(426) 



(427) 



k>0 



where N l and the antisymmetric M y define constant tensors. 

• Let us now construct the rank-2 Killing tensors which solve the Killing equation, 

djC kl (x, tf) + d t C jk (x, $ + d k C lj (x, = . (428) 

Considering the expansion in terms of Grassmann degrees of freedom [Berezin] of the Killing 
tensor C J (x, ip) , we get the coefficients tensors Q{... afc which are constructed as sym- 
metrized products [Gibbons 1987] of Yano-type Killing tensors, C l Y {x) , associated with the 



rank-1 Killing tensors C l (x) obtained by (426) 

1 



(429) 



It is worth noting that the Killing tensor defined in (429) is symmetric in its bosonic indices 
and anti-symmetric in the fermionic indices. Thus, we obtain 

C«(x, J) = J2 (M ( l n Wl n x l x m + M^N^x 1 



k>0 



+N^M j l n x m + N$Nfi n 
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(430) 



i\) ax ■ ■ ■ tp ak 



where M % \ , M % \ , N J k and iVj? are skew-symmetric constants tensors. Then one can verify 
with direct calculations that (427) and (430) satisfy Killing equations. 



Having constructed the generic Killing tensors (427) and (430) generating constants of 



the motion, we can now describe the supersymmetries of the Pauli-like Hamiltonian (412) 



To start, we search for momentum-independent invariants, i.e. which are not derived from 
a Killing tensor, C % = C l i = ■ ■ ■ = . In that event, the system of equations (423) reduces 
to the two constraints, 



dQ c (x,ip) 

# or Fam = ' o(0) 

diQ c (xJ)=0. o(l) 



(431) 



For (7 = 0, which means no spin-gauge field coupling, it is straightforward to see that 
the spin vector, in particular, and all arbitrary functions / (ip ) which depend only on the 
Grassmann variables are conserved along the motion. 

For nonvanishing gyromagnetic ratio g, only the "chiral" charge 



(432) 



remains conserved. The "chiral" charge Q c can be considered as the projection of the 
internal angular momentum, S, onto the internal trajectory ip{ T ) ■ Thus, Q c can be viewed 
as the internal analogue of the projection of the angular momentum, in bosonic sector, onto 
the classical trajectory x(t) . 

Let us now construct superinvariants linear in the covariant momentum. — ■ ■ ■ 



such that (423) becomes 



leg 



e 9 im dC (^^) 



e 9 im dC j (x,ip) 



d,C(x^) = eF jk C k (x^)-^ m 
^ djC k (x , if) + d k C j (x, tf) = . 



0, o(0) 

o(l) 
o(2) 







(433) 



Choosing the non-vanishing term N 3 a = 5 J a , in the general rank-1 Killing tensor (427), 



provides us with the rank-1 Killing tensor generating the supersymmetry transformation, 



ci(x,ij) = 5ir 



(434) 
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By substitution of this Grassmann- valued Killing tensor into the first-order equation of (433) 
we get 



(435) 



Consequently, a solution C(x,ip) = of (435) is only obtained for a fermion with ordinary 
gyromagnetic ratio 

9 = 2. (436) 

Thus we obtain, for V(r) = 0, the Grassmann- odd supercharge generating the M = 1 
supersymmetry of the spin-monopole field system, 



Q = ^-I1, {Q, Q} = -2iH 2 . 



(437) 



For nonvanishing potential, V(r) 7^ , the zeroth-order consistency condition of (433) is 
expressed as 13 



V r = 



(438) 



Consequently, adding any spherically symmetric potential V(r) breaks the supersymmetry 
generated by the Killing tensor C J = 5 3 a ip a : M = 1 SUSY requires an ordinary gyromagnetic 
factor, and no additional radial potential is allowed |Spector| . 



Another Killing tensor deduced from (427) is obtained by considering the particular case 



with the non-null tensor Ni 



a\a,2 



This leads to the rank-1 Killing tensor, 



(439) 



In this case, the first-order constraint of (433) is solved byC(x, if) = 0, provided the gy 



romagnetic ratio takes the value g = 2 . For vanishing potential, it is straightforward to 
verify the zeroth-order consistency constraint and therefore to obtain the Grassmann- even 
supercharge, 

Qi = £-n, (440) 
defining the "helicity" of the spinning particle. As expected, the consistency condition of 



superinvariance under (440) is again violated for V(r) 7^ , breaking the supersymmetry of 



the Hamiltonian , in (437). 



We use the identity S k G j d B k = V>V m G^djF lm = . 
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Let us now consider the rank-1 Killing vector, 



C j (x,$) = (Sxx) j 



(441) 



obtained by putting M % 1 = (i/2)e 3 €k ai a 2 into the generic rank-1 Killing tensor (427) 



The first-order constraint is satisfied with C(x,ip) = 0, provided the particle carries gyro- 
magnetic ratio g = 2. Thus, we obtain the supercharge, 



Q 2 = {xxU)-S. 



(442) 



which, just like those in (437) and (440) only appears when the potential is absent, V = 0. 



We consider the SUSY given when M l3 a = e a u so that the Killing tensor (427) reduces to 



C'(x^) 



(443) 



-> Q ifj ■ X 

The first-order constraint of (433) is solved with C(x, ift) = — (g — 2) . The zeroth-order 



consistency condition is, in this case, identically satisfied for an arbitrary radial potential. 
We have thus constructed the Grassmann-odd supercharge, 



Q3 = (ifxfl).^+!(0-2)^ 
2 r 



(444) 



which is still conserved for a particle carrying an arbitrary gyromagnetic ratio g . Note, 
that this supercharge generalizes the one obtained in the restricted case with g = 2 
DeJonghe 1995] . See also |Horvathy 20001 . 



Now we turn to invariants which are quadratic in the covariant momentum. For this, we 
solve the reduced series of constraints, 

e 9.,.m dC r. 



C%V + ^ ^ m C 3 d 3 F lm - ^ m j—F am = 0, 
4 2 oip a 

a-C = C jk d k V + eF 3k C k + ^VV™ C jk d k F lm - ^ m<K ' 



dijj a 



en dC jh 

d 3 C k + d k C* = e (F ]m C mk + F km C m >) - y #™^r F * 



o(0) 

o(l) 
o(2) 
o(3) 



(445) 



djC km + d m C jk + d k C mj = . 
We first observe that C lJ (x,ip) = 5^ is a constant Killing tensor. Solving the second- and 



the first-order constraints of (445), we obtain 



C j (x,^) = and C(x,ip) = V(r) 
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eg 



S-B 



(446) 



respectively. The zeroth-order consistency condition is identically satisfied so we obtain the 
conserved energy of the spinning particle, 



f = ^if 2 - e 4s-B + V(r) 



2~ 2 < 447) 

Next, we introduce the nonvanishing constants tensors, M % i k = e^ k , N lJ a = , into 



(430) in order to derive the rank-2 Killing tensor with the property, 

C jk (x, $) = 2 5 jk {x ■ if) - x J if) k - x k if) j . 



(448) 



Injecting the Killing tensor (448) into (445), we satisfy the second-order constraints with 

if) x x 



C(x,iP) 



Q 



9 (2-<7) 
2 r 



(449) 



To deduce the integrability condition of (445), we require, in the first-order constraint, the 
vanishing of the commutator, 



[d h d j }C(x)=0 



A(V(r)-(2-g) 



2 <f 



0. 



(450) 



Then the Laplace equation (450) provides us with the most general form of the potential 



admitting a Grassmann-odd charge quadratic in the velocity, namely with 

,2 
~2 



V(r) = (2-g) 2 ^- 2 + - + p. 

Sr z r 



Consequently, we solve the first-order constraint with 



fa 



C(x, if)) = y egS ■ B ) x • if) . 



(451) 



(452) 



so that the zeroth-order consistency constraint is identically satisfied. Collecting our results 
leads to the Grassmann-odd conserved charge quadratic in the velocity Ngome 03/2010] , 



Q 4 = (fix (fxfi)) -^+|(2- g)^^ (--egS-B^jx-ijJ. (453) 

Let us underline that the conserved charge Q 4 which is not a square root of the Hamiltonian 
Tig remains conserved without restriction on the gyromagnetic factor, g . We can also remark 
that for g = , this charge coincides with the scalar product of the separately conserved 
Runge-Lenz vector 14 [Mcintosh 1970} Zwanziger 1968] by the Grassmann-odd vector: 



Q 4 



9=0 



K 



s=0 • 



ifj. 



(454) 



14 



The case of spinning particle with null gyromagnetic ratio, g = , coincides with a spinless particle. 
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The supercharges Q and Qj with j = 0, • • • , 3 , previously determined, form together, 
for ordinary gyromagnetic ratio, the classical superalgebra, 

{Qo, Qo} = {Qo, Qi} = {Q, Qi} = {Qi, Qx} = {Q 2 , Q 2 } = o, 
{Qo,Q} = iQu {Qo, Q 2 } = {Q 2 , Q 3 } = 0, 

{Qo, Q 3 }=*Q 2 , {Q, Q} = -2m 2 , (455) 
{Q, Q 2 } = {Qi, Q 3 } = Q 4 , 

{Q, Q 3 } = 2iQi, {Q 1; Q 2 } =zQ 3 Q, {Q 3 , Q 3 } = i (2Q 2 - Q 5 ) , 

where Q5 is a bosonic supercharge that we will construct below |465 . From (455) it follows 
that the linear combination Qy = Q 3 — 2Q has the special property that its bracket with 
the standard supercharge Q vanishes: 



{Qy, Q} = 0. 



(456) 



Indeed, Qy is precisely the Killing- Yano supercharge constructed by De Jonghe, Macfarlane, 
Peeters and van Holten |De Jonghe 19951 - 



Let us now investigate the bosonic symmetries of the Pauli-like Hamiltonian (412). We 



use the generic Killing tensors previously constructed [ cf. (427) and (430)] to derive the 
associated bosonic constants of the motion. 

Firstly, we describe the rotationally invariance of the system by solving the reduced series 



of constraints (433). For this, we consider the Killing vector provided by the replacement, 



= —e t:l k n k into (427). Thus for any unit vector n, we obtain the generator of space 



rotations around n. 



C(x, ip) = n x x . 



Inserting the previous Killing vector in the first-order equation of (433) yields 



C{x,ij}) = c{4>) - q 



n ■ x 



(457) 



(45£ 



Moreover the zeroth-order consistency condition of (433) requires for arbitrary radial poten- 
tial, 

c(tf) = S-n. (459) 
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Collecting our results provides us with the total angular momentum, which is plainly con- 
served for arbitrary gyromagnetic ratio, 

j = L + S = xxli-q- + S. (460) 

r 

In addition to the typical monopole term, the conserved angular momentum also involves 
the spin vector, S . It generates an o(3) rotations bosonic symmetry algebra, 

{j\J j } = e ijk J k . (461) 

In the particular case of vanishing gyromagnetic factor g = 0, the usual monopole angular 
momentum L and the internal spin angular momentum S are separately conserved involving 
an 

o(3) 

rotations © o{3) spin (462) 

symmetry algebra . 



We turn into invariants which are quadratic in the covariant momenta. Then, we have 



to solve the series of constraints (445). We first observe that for M^ mk — M^ mk — e^ mk , the 



Killing tensor (430) reduces to the rank-2 Killing-Stackel tensor, 



C ij (x,$) = 2S ij x 2 -2x i x j . 



(463) 



Inserting (463) into the second- and in the first-order constraints of (445), we get for any 



gyromagnetic factor and for any arbitrary radial potential, 



C(x,ip) = and C(x,ip) = —gq 



x ■ S 



(464) 



Hence, we obtain the Casimir 



Q 5 = J 2 -q 2 + (g-2) J-S-gQ 2 . 



(465) 



The bosonic supercharge Q 5 is, as expected, the square of the total angular momentum, 
augmented with another, separately conserved charge Ngome 03/2010] , 



(g-2)J-S-gQ 2 . 



(466) 



Indeed, for g = 0, (466) directly implies that the product, J ■ S , and hence the spin 



vector, 5* , are separately conserved. 
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• For g = 2 , we recover the conservation of the supercharge Q2 [ cf. (442) . 

• For the anomalous gyromagnetic ratio g = 4 , we obtain that J ■ S — 2Q 2 is a constant of 
the motion. 



Now we are interested in the hidden symmetry generated by a conserved Laplace-Runge- 



Lenz-type vector. Then, we introduce into the algorithm (445) the generator 



C tj (x, ip) = 2 6 lJ ft ■ x - n V - n j x l 



(467) 



easily obtained by choosing the non-vanishing, iV* J = e im % m and M % i m = e^™ 1 , into the 



generic rank-2 Killing tensor (430). Inserting (467) into the second-order constraint of 



(445), we get 



C(xJ) = q r ^ + C(j). 



(468) 



We solve the first-order constraint of (445) by expanding C(x, ip) in terms of Grassmann 
variables [Berezin] . 

C(x, if) = C(x) + Car-* W 11 • • • r k ■ (469) 



k>l 



Consequently, the first- and the zeroth-order equations of (445) can be classified order-by 



order in Grassmann-odd variables. Thus, inserting (468) in the first-order equation, and 
requiring again the vanishing of the commutator, 

,2 



[di, d 3 }C(x)=0 



A V(r) 



2r< 







(470) 



we deduce the most general radial potential admitting a conserved Laplace-Runge-Lenz 
vector in the fermion-monopole interaction, namely 



V{r) 



2r 2 r 



+ 7, /i ,7 G R. 



(471) 



Investigating the first term on the right-hand side of (469), we obtain 

C [x) = fi . 



(472) 



Introducing (468) and (471) into the first-order constraint of (445), on one hand, provides 
us with 



(473) 
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and on the other hand with 



k>l 



eg 



S-B)(n-x)- 



gq 



1 - 



g\n-S 



(474) 



with g(g - 4) = . 



Let us precise that the zeroth-order consistency condition of (445) is only satisfied for 



C(i/>) 



-S ■ n . 



Q 



(475) 



Collecting our results, (467), (468), (471) and (474), we obtain a conserved Runge-Lenz 
vector if and only if 

g = or g = 4 ; (476) 

we get namely 



x 



n x j + fji- + l - - )s x n- 



eg 



S-B)x 



gq 

2 



1 - 



9\ S u 



2/ r 



q 



(477) 



Note that the spin angular momentum which generates the extra "spin" symmetry for 
vanishing gyromagnetic ratio is no more separately conserved for g = 4. Then, an interesting 
question is to know if the extra "spin" symmetry of g = is still present for the anomalous 
superpartner g = 4 in some "hidden" way. 

Let us consider the "spin" transformation generated by the rank-2 Killing tensor with 
the property, 

~* 9 I c""„fc i ck^m\ (478) 



C mk [x , V) = 25 mk {S-n)- u - (S m n k + S k n m ) . 



The rank-2 Killing tensor (478) which can be separated as C mk = C™ k + C™ k is obtained 
by putting 



iki 



,_ f jk f m 
2 m a\a,2 ' 



jkl m 
c "m c ai a 2 i 



(479) 



into the general rank-2 Killing tensor ( |430 ). Inserting (478) into the second-order constraint 



of (14451) leads to 



C(x , tf) = ( S X U > + (7(V>) and ^(^-4) = 0^ 



(480) 



We use the potential ( 471| to solve the first-order equation of (445J), 

S ■ n + c(-0) . 

p(p-4) =0. 



C(^^=(2V(r)-^ 



4r 



2g 



(481) 



and 
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The zeroth-order consistency condition is satisfied with 

c^) = ~ 9 4^S-n, (482) 
o q 

so that collecting our results provides us with the conserved "spin" vector, 

^=(n 2 + (2-^)^))5-|(n.^n + f(^ + ^)5xn 



"T^ -7 ^ with 9(9-^) =o. 



(483) 



In conclusion, the additional o(3) spin "spin" symmetry is recovered in the same particular 



cases of anomalous gyromagnetic ratios and 4 [cf. (476) 

• For g = 0, in particular, 

n = 2SS. (484) 

• For g = 4, we find an expression equivalent to that of D'Hoker and Vinet 



D'Hoker 09/1985] , namely 

(l 4 = (^n 2 -2K(r))5-2(n-5)n + 2(^ + ^)5xIl-4^-7^ S. (485) 

Note that this extra symmetry is generated by a Killing tensor, rather than a Killing vector, 
as for "ordinary" angular momentum. Thus, for sufficiently low energy, the motions are 
bounded and the conserved vectors J , K g and Q g generate an 

o(4) © o(3) spm (486) 

bosonic symmetry algebra. 

So far we have seen that, for a spinning particle with a single Grassmann variable, SUSY 
and dynamical symmetry are inconsistent, since they require different values for the g- 
factor. Now, adapting the idea of D'Hoker and Vinet to our framework, we show that the 
two contradictory conditions can be conciliated by doubling the odd degrees of freedom. 
The systems with g = and g = 4 will then become superpartners inside a unified M = 2 
SUSY system jFeher 1988] . 

We consider, hence, a charged spin-^ particle moving in a flat manifold Ai D+2d , inter- 
acting with a static magnetic field B . The fermionic degrees of freedom are now carried 
by a 2<i-dimensional internal space |Belluccil IKochan| IGonzalest Avery 2008 . This is to be 
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compared with the d-dimensional internal space sufficient to describe the M = 1 SUSY 
of the monopole. In terms of Grassmann-odd variables ipi % , the local coordinates of the 
fermionic extension Ai 2d read [if)f, ip^) with a,b — 1, • • • , d. The system is still described 



by the Pauli-like Hamiltonian (412). Choosing d = 3 , we consider the fermion £ Q which is 



a two-component spinor, £ a = ( 1 ] , and whose conjugate is £ Q [Salomonsonj . Thus, we 

W 

have a representation of the spin angular momentum, 

S k = \l a v k Jtp with a,/? = 1,2, (487) 

where the a k a ^ with k = 1,2,3 define the standard Pauli matrices. Defining the covariant 
Poisson-Dirac brackets as 



+ *(-i) a/ ( 



df dh + df dh 



(46 



(489) 



we deduce the non- vanishing fundamental brackets, 

{x 4 ,^} = 5j, {n,,^} =ee i3k B\ = -iSf, 

{S k ,S 1 } = e kl m S m , {S k } ^} = - Z -ea k /, {S k ,^} = % -aYi u . 

We also introduce an auxiliary scalar field, $(r) , satisfying the "self- duality" or " Bogo- 
molny" relation 15 , 

{lT\$(r)} = ±eB k . (490) 

This auxiliary scalar field also defines a square root of the external potential of the system 
so that 

l$ 2 (r) = V(r) . (491) 



As an illustration we obtain the potential 16 defined in (471) by considering the auxiliary 
field 

d>(r) = ±{l + f)- (492) 



See |Feher 1988] to justify terminology. 

The constant is 7 = — - . 

' 2q 2 



lo 
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In order to investigate the M = 2 supersymmetry of the Pauli-like Hamiltonian (412), 
we outline the algorithm developed we use to construct supercharges linear in the gauge 
covariant momentum, 



3 _|_ r)k I ^kv 



dC dC 



11 de 



f ev^CidjB^O, o(0) 



3 m C = ee mjk B«C> + &B k ( fa*"— - — 



4 



d j C h (x,Z,£) + d k &{x i Z,g) = 0. 
Let us first consider the Killing spinor, 



o(l) (493) 
o(2) 



(494) 



Inserting this Killing spinor into the first-order equation of (493) provides us with 

d m Cp = --e B m £ p and g = 4, 



(495) 



which can be solved using the self-duality relation (490). Thus, we get 



(496) 



provided the anomalous gyromagnetic factor is g — 4. The zeroth-order constraint of (493) 
is identically satisfied, so that collecting our results provides us with the supercharge, 



Q^ = ^n,<^±^(r)^ 



(497) 



To obtain the supercharge conjugate to (497|), we consider the conjugate Killing spinor, 



2 q a 



(498) 



In the case of anomalous value of the gyromagnetic ratio g = 4, the first-order equation of 



(493) is solved by using the Bogomolny equation (490). This leads to the conjugate 



(499) 



The zeroth-order consistency constraint is still satisfied, so we obtain the odd-supercharge, 

Q^ = lCa k JU kT ^(r)^. (500) 
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The supercharges Qp and Q 13 are, both, square roots of the Pauli-like Hamiltonian "H 4 , 

{Q p ,Qp} = -iHtl, (501) 

and therefore generate the M = 2 supersymmetry of the spin-monopole field system. It is 
worth noting that defining the rescaled, 



1 



Ha 



and Ur 



1 



it is straightforward to get, 



7^0 — U P 7^4 Ur , 



(502) 



(503) 



which make manifest the fact that the two anomalous cases g = and g = 4 can be viewed 
as superpartners 17 , see |Feher 1988] . Moreover, in our enlarged system, the following bosonic 
charges 



x 



J = x xll - q- + S 



K = II x J + fi - - S x II - 2e (S ■ B) x + 2q- + -S , 
r V / r q 

n = Q? 3° Qa = \ ($ 2 (r) - n 2 ) s + (n • s) ff t $(0 s x n, 



(504) 



remain conserved such that they form, together with the supercharges Qs and , the 



classical symmetry superalgebra D'Hoker 09/1985 IFeher 1988] 



{G?,Qp} = -iHtl, {Q",&} = {Q,,Q,} = 0, {Q^,J k } = - A Q a a k J 



1 (t n*„iP 
4 q 



4 q 



{K\K j } = e 



ijk 



2U A J k + 2Q* 



{n\Q?} = e ijk u±n k , 



{ T , A J } = e ijk A k with A 1 = J 1 , K l , SI 1 . 

We have shown, in this section, that the Runge-Lenz-type dynamical symmetry and the 
additional extra "spin" symmetry both require instead an anomalous gyromagnetic ratio, 



g = or g — 4 . 



(505) 



17 With The scalar = 2 . 
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These particular values of the g-factor come from the effective coupling of the form Fij =F 



6ijkD k §, which add or cancel for self-dual fields |Feher 1 988] . 



e ijk D k $ . 



(506) 



Moreover, the super- and the bosonic symmetry can be combined in this enlarged fermionic 
space and provides us with an TV = 2 SUSY, as proposed by D'Hoker and Vinet 
D'Hoker 09/1985] . See also jFeher 19881 IFeher 1989} |Feher 02/1989] IBlooreJ. 



At last, let us remark that confining the spinning particle onto a sphere of fixed radius 



p implies the set of constraints DeJonghe 1995 



x 



P 



x ■ if) — and x ■ H — . 



(507) 



This freezes the radial potential to a constant, and we recover the TV = 1 SUSY described 
by the supercharges Q , Q\ and Q2 for ordinary gyromagnetic factor g = 2 . 
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B. N = 2 SUSY in the plane 



The planar system consisting of a spinning particle interacting with a static magnetic 
field in the plane exhibits more symmetries as its higher- dimensional counterpart. Indeed, 
the M = 2 supersymmetry, here, is realized without doubling the Grassmann-variable of the 



internal space as it was the case in three-dimensional space system, see section V A Such 
an "exotic" supersymmetry, which is realized in two different ways, is only possible in two 
spatial dimensions [Duval 1993} IDuval 19951 IDuval 2008] . This is one more indication of the 
particular status of two-dimensional physics. 

To see this, we investigate the two dimensional model given by the Pauli-like Hamiltonian, 

H = \ ^ ~ ~ SB + V ^ ' (508) 
where the magnetic field simplifies into 18 

Fij = e i3 B = diAj - djAi , (509) 

and the spin tensor is actually a scalar 



S = - % -e i pl> i V- (510) 



The fundamental brackets remain the same as in (417), and the spatial and the internal 



motions of the particle are governed by the following equations, 

go 

x k = —Sd k B + eBe kj x J + d k V , 

(511) 

^ = ^8^, S = 0. 

Observe the conservation of the spin S along the particle motion and let us recall that all 
quantities quadratic in the Grassmann variables are proportional to S. 



We search for dynamical quantities which are constants of the motion, for the planar 



IS 



We dropped the irrelevant third z -direction. 
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system, by solving the series of constraints: 



o(0) 
o(l) 



dipi dipi 

&H. dC- 

diCj + djCi = e (F ik C k j — C ik F k j) + i^—^-^- + Cij^d^H , o(2) 

diC jk + djCki + d k Cij = CijkidiH + (terms linear in Ci mn ) . o(3) 



(512) 



Using the equality 



.an 



eg 



2 J 



eg 
' 2 



the zeroth-order constraint in (512) becomes 

eg „ dC 
~2 



eg 



(', [diV - -j-SdiB 



complemented by the first-order equation of (512) 



g 



ac, 



diC = eB ( e l3 C 3 + - e jk ^ 



Similarly the second and higher-order equations take the form 

= eB (^EikCkj + £jkCki + ^ £ i k ^i~Q^jf S j (pkV — ^- SdkB 



For radial functions V(r) and B(r) 



hence 



diV = —V. 
r 



eg 



diB = X -± B ' 
r 



d s V - ■;• Sd,B ) C,., = ^ ( V* - e 4 SB') C L .,. 



(513) 



(514) 



)+C ij (d j V- e ^ SdjB). (515) 



(516) 



(517) 



(518) 



r v 2 

Let us now consider some specific cases. To this, we introduce the universal generalized 
Killing vectors, 

Ci = EijX 3 , Ipi, Eijlp 3 } , (519) 

where 7, denotes a constant vector. 

• A constant Killing vector 7^ gives a constant of the motion only if we can find solutions 
for the equations 



dC ( 2 

diC = eBei^j, Bejiipi— = 7; ^— diV - SdiB ) . 



(520) 
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Now for a Grassmann-even function 



C = c + c 2 S . 



(521) 



the left-hand side of the second equation in (520) vanishes, therefore we must require B and 
V to be constant. This leads to the solution 



C = -eBe ij ^ l x j 



V = const, B = const . 



(522) 



The corresponding constant of the motion, ( , is identified with the "magnetic translations" 
Hughes 1986] , 

(523) 



( = YPi with P, t = Ui - eBe^x 3 



Next we consider the linear Killing vector Cj = £«x J , with all higher-order coefficients 



Cij... = 0. Again for Grassmann-even C the left-hand side of equation (514) vanishes, and 
we get the condition 

(524) 



which is automatically satisfied for radial functions B(r) and V(r). Therefore we only have 



to solve the equation (515): 



diC = -eBxi = — — (rB). 
We infer that C{r) is a radial function, with 

C = -erB . 



(525) 



(526) 



Therefore C is given by the magnetic flux through the disk D r centered at the origin with 
radius r: 

72^- e * ( 527 ) 



C = -— B(r)d 2 x = -—^ B (r). 

27T J D 27T 



We then find the constant of the motion representing the angular momentum 
Ngome 03/20101 , 

(528) 



associated with the o(2) 



rotations 



L = Eijxjlj + — $ B (r), 
Ztt 

symmetry group. 



Ill 



• There are two Grassmann-odd Killing vectors, the first one being C, = ipi- With this 
Ansatz, we get for the scalar contribution to the constant of the motion the constraints 



eg 



dC 



eB 



2 B Eijipj ^— = ipidiV and d { C = — (2 - g) e^ipj 



(529) 



It follows that either g = 2 and [C, V) are constant, in which case one may take C = V = 0, 
or g 7^ 2 and C is of the form 



C = EijK^r)^ with dy = - e 4~ BK, 



diKj= {2 ° )eB 5 ir (530) 



This is possible only if B is constant and 
_ eB(2-g) = 



K(rJ = e B r 



It follows that we have a conserved supercharge of the form, 



egn 
Air 



(531) 



Q = i) 1 (Ik - KEijX j ) 



(532) 



The bracket algebra of this supercharge takes the form 



i{Q,Q} = 2U + (2-g)eBJ, J = L + S. 



(533) 



Of course, as S and L are separately conserved, J is a constant of the motion as well. It is 



now easy to see that for ordinary g -factor the supercharge Q in (532) is a square root of 
the Hamiltonian. 

• Let us remark that for the anomalous gyromagnetic ratio g — 1 , we construct the 
conserved conformal supercharge, 



S = x-if>-tQ, 



(534) 



obtained by using the internal equation of the motion in (511) 



Finally we consider the dual Grassmann-odd Killing vector Cj = Sijijjj. Then the 



constraints (514) and (515) become 



B — — = diV, diC = ip h 

2 atpi 2 



implying that 



C = N i {x)1> i and ^-BN l = d t V, d i N j = (g ^ eB 



(535) 



(536) 
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As before, the magnetic field B must be constant and the potential is identical to (531) 



N^-KXi, V = -^® B {r) = 9 -^—^e 2 B 2 r 2 . (537) 
Thus, we find the dual conserved supercharge 19 |Horvathy 2 005 1, 

Q* = Eijlpi (Uj - KE jk X k ) = Xpi {SijU-j + KXi) , (538) 

corresponding, in the case of ordinary gyromagnetic ratio, to the "twisted" supercharge 
used by Jackiw |Jackiw 1984j to describe the Landau states in a constant magnetic field. 
Moreover Q* satisfies the bracket relations 

i{Q*,Q*} = 2H + (2-g)eBJ, i{Q,Q*} = 0. (539) 

Thus the harmonic potential ( 531[ ) with constant magnetic field B allows a classical M = 2 



supersymmetry with supercharges (Q, Q*) , whilst the special conditions g — 2 and V = 
allows for M = 2 supersymmetry for any B(r). 

• As a consequence of the conservation of the "twisted" supercharge, we construct for 
g = 1 , the associate conserved conformal supercharge 

S* = x x ijJ+tQ*. (540) 

Thus, for the non-ordinary gyromagnetic ratio g = 1 , the supercharges Q , Q* , S and 
S* extend the o(2, 1) algebra into an osp(l,l) superalgebra and satisfy the commutation 
relations, 

{Q, Q} = {Q\Q*} = -i{2H + eBJ) , {S\S*} = -it 2 (2n + eBJ) , 
{Q*,S*} = -{Q,S} = -it^H + eBj) +ix -n, {Q,Q*} =0, 

(541) 

{Q*,S} = {Q,S*} = -i(L + 2S) , {S,S*} =2it(L + 2S) , 

{5,5} = -it 2 (2H + eBJ) + lit x ■ fl - ir 2 . 
The van Holten recipe is therefore relevant to study planar fermions in an arbitrary planar 
magnetic field, i.e. one perpendicular to the plane. As an illustration, we have shown, for 



ordinary gyromagnetic factor, that in addition to the usual supercharge (532) generating 
the supersymmetry, the system also admits another square root of the Pauli Hamiltonian 
% |Horvathy 20 05 . This happens due to the existence of a dual Killing tensor generating 
the "twisted" supercharge. 



19 



The cross product of two planar vectors, a x b — Eija l W , again defines a scalar. 
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VI. NON-COMMUTATIVE MODELS 



A non- commutative oscillator with no kinetic term but with a certain momentum- 
dependent potential is constructed. The classical trajectories followed by a non-commutative 
particle in this oscillator field lie on (arcs of) ellipses. 

A. Non- commutative oscillator with Kepler-type dynamical symmetry 

In recent years, a remarkable non-commutative model was derived in the context of solid 
state physics by Chang and Niu Chang 1995| . They stated that the semiclassical analysis 



of a Bloch electron in a three-dimensional crystal lattice reveals an extra "Berry phase" 
term, , which can take a monopole-like form in the band structure. The study of the 
wave-packet dynamics of this Bloch electron, under perturbations slowly varying in space 
and in time, leads to the equations of the motion in the mth band [in units h — 1], 

k = -eE-eS x B(x) , x*= d£m {^ -kxQ(k). (542) 

dk 

Here £ m (k) , x and k denote the Bloch electron's band energy, the intracell position and 



the quasi-momentum, respectively. Note that in the right hand side equation of (542), the 
electron velocity gains an anomalous velocity term, k x Q(k) , which is the mechanical 
counterpart of the anomalous current. 



In a magnetic field-free theory [with B = 0], the equations (542) can also be deduced 
using the symplectic closed two-form, 

Q = dpi A dxi + ^eijk® l dpj A dpk , (543) 

where the "extra" term induced by the Berry phase yields the position coordinates non- 
commutative |Chang 1995[ INiuj , 



{xi,Xj} = e ijk Q k = Qij, {%i,Pj} = Sij, {Pi,Pj} = . (544) 



Applying the Jacobi identities to the coordinates, we get 

d& m 

U \Pii V^j'j %k / J cyclic £jkm g^. ' 

QQij QQjk QQki 
U \ \Xj } Xk j j cyclic n ~n 

dp k dpi dpj 
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(545) 



Then, the vector field has the property 20 to be only momentum-dependent [Berard 20 04] , 

e i = e i (p), (546) 

and also requires the consistency condition 

V p --G(p)=0, (547) 

which can be interpreted as a field Maxwell equation in the dual momentum space. Choosing, 
for example, the non-commutative vector aligned in the third direction, 



@i = 96 i3 , 9 = const , 



(548) 



the 3-dimensional theory reduces to the planar mechanics based on "exotic" Galilean sym- 
metry jLukierski 19971 IDuval 20001 IDuval 20011 IChaichian 20"0T1 INair 20011 IScholtz 20051 



IScholtz 2 009. H orvathy 2010| . As an application of ((5481, some interesting results, includ- 



ing perihelion point precession of the planetary orbit, can be derived [Romero 2003] when 
taking into account the Kepler potential, 



V(r) 



oc r 



(549) 



Other applications of (548) concern, for example, the Quantum Hall Effect [Dunne 19901 
IDuval 2000| |Horvathy 2"002 . 



Such a choice only allows for axial symmetry, though. In our theory, however, we re- 
store the full rotational symmetry by choosing instead G to be a "monopole in p-space" 
[Berard 2004] . 



G, = 9 



TEH 
P 3 



9 = const 



(550) 



where p — \p\ ■ Indeed, away from the origin, the dual monopole (550) is the only spherically 



symmetric possibility consistent with the Jacobi identities 21 . Let us mention that the p- 



monopole form in (550) has already been observed experimentally by Fang et al. in the 



context of anomalous Hall effect in the metallic ferromagnet SrRu03 |Fang 2 003 1. 



As expected, (550) corresponds to extra, "monopole" term in the symplectic structure 



(543 ) which is in fact that of a mass-zero spin-# coadjoint orbit of the Poincare group. The or- 



bit is indeed that of the o(4, 2) conformal group [Penrose 1 9 72 \ [Penrose 19771 ICordani 19 90] . 



20 For a more general theory which also includes magnetic fields, see, e.g., Chang 1995 Niu. Duval 2000 . 
For simplicity, the mass has been chosen unity. 



21 



See the equivalent demonstration in real x-space in section III A 
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We can now study the 3D mechanics with non-commutativity (550), augmented with the 
Hamiltonian, 

n = ^ + V(x,p), (551) 

where we allowed that the potential may also depend on the momentum variable, p 22 . 
The equations of motion of the system read 

Xi=Pi + ^- + 9eij]M^-, Pi = -jr-> ( 552 ) 
dpi p A oxj axi 



where, in the first relation, the "anomalous velocity terms" is due to our assumptions (550). 

We are particularly interested in finding conserved quantities. This task is conveniently 
achieved by using van Holten's covariant framework |van Holten 20 07J. which amounts to 
searching for an expansion into integer powers of the momentum, 

Q = C (x) + Ci(x)pi + ^CijfflpiPj + ^CijkfflpiPjPk + ... . (553) 

Requiring Q to Poisson-commute with the Hamiltonian yields an infinite series of con- 
straints. However, the expansion can be truncated at a finite order n , provided to satisfy 
the Killing equation, -D^C^...* \ = , when we can set Cj 1 ...j n+1 ... = . 

Let us assume that the potential has the form V(|a?|, \p\) , and try to find the conserved 
angular momentum, associated with the Killing vector C = n x x , which represents space 
rotations around n. An easy calculation shows that the procedure fails to work, however, ow- 
ing to the p-monopole term. We propose, therefore, to work instead in a "dual" framework 
Ngome 06/2010] , i.e. in momentum space, and search for conserved quantities expanded 



rather into powers of the position, 



Q = C (p) + Ci{p)Xi + ^Cij(p)xiXj + ^Cijkip)XiXjX k .... (554) 



22 Note that momentum-dependent potentials are frequently used in nuclear physics and correspond to 
non-local interactions. 
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Then, the covariant van Holten algorithm, presented in section II B is replaced by 

= o(0) 
= o(l) 

0^ (e ijm C lk + e lkm C i3 ) - + + C ijk (pi + =0 o(2) 



r or \ p A opj J \ opt 

1 dV ( n p, 
r Or 

ldV ( pm, r r r s fdCij dC jk dC ki A 

r Or \ p 6 v Opk opi Opj J 

dV 



Cujk I pi + 



dpi 



o(3) 



where r = \x\. The expansion (554) can again be truncated at a finite order n , provided the 



higher order constraint of the previous series of constraints transforms into a dual Killing 
equation, 



d(p H C Pi2 ..^ n) - 0. 



(555) 



Then, for linear conserved quantities, Q = C (p) + Cj(p)xj , we can set CV,- = C^k = . . . 0. 
The dual algorithm therefore reduces to 



a(p i+ -yo, o(o) 

< eei ^ Ci -wr 0, 0(1) 

dC k dC j 

0(2) 

Introducing the dual Killing vector 

C = n x p 

into the previous algorithm provides us with 



C = On ■ p . 



P 



P 
P 



Thus, we obtain the conserved angular momentum, 



J = L 



x x p — 6p . 



(556) 



(557) 



(558) 



which is what one would expect, due to the "monopole in p -space", whereas the non- 
commutative parameter, 9, behaves as the "monopole charge" [Cortes 1 996]. 
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The next step is to inquire about second order conserved quantities. Then, the series of 
constraints which has to be solve read 



o(0) 



- B ^(e^- d f)+c ii L + B f)= o,o(i) 

r or \ p 6 opj J \ dpi J 



( e ij m Cik + tikmPij) — (-75 V — , o(2) 



dC k dd 



(559) 



p 

dCu dC 



dp, dp k 



+ 



+ 



dC, 



hi 



0. 



0(3) 



k dp k dpi dpj 

Remark that usually the Runge-Lenz vector is generated by the rank- 2 Killing tensor CV, = 
IbijU ■ x — riiXj — rijXi where n is some fixed unit vector |van Holten 2007] . Not surprisingly, 
the original procedure fails once again. The dual procedure works, though. The dual two- 
tensor 

Cy = Ibijn ■ p- mpj - rijPi , (560) 



verifies the dual Killing equation of order 3 in (559). Then the order-2 equation yields 



C = 9 



n x p 
p 



(561) 



Inserting into the first-order constraint of (559) and assuming d r V 7^ 0, the constraint is 
satisfied with 

C = an-p (562) 

a being an arbitrary constant, provided the momentum-dependent potential and the Hamil- 
tonian take the form 



x 2 p 2 9 2 a . „, x 2 9 2 a 

V = — H + - and H = 1 ~ + - , 

2 2 2p 2 p 2 2p 2 p ' 

respectively. Then the dual algorithm provides us with the Runge-Lenz-type vector 



(563) 



K = x x J — ap . 



(564) 



Its conservation can also be checked by a direct calculation, using the equations of the 
motion, 



x = 9 



x x p 

p3 



9 2 a 

"4 + T J P: 



p = —x 



(565) 



where the anomalous velocity term in the first relation is transversal. 
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Note that the (—p 2 /2) term in the potential cancels the usual kinetic term, and our sys- 
tem describes a non-relativistic, non- commutative particle with no mass term in an oscillator 
field, plus some momentum- dependent interaction. 

Writing the Hamiltonian as 

x 2 9 2 fl a\ 2 a 2 ,„ N 

u = Y + ^{ P + ¥) ~W 2 (566) 
a 2 e 2 

shows, moreover, that % > — with equality only attained when p = , which plainly 

20 1 a 

requires a < 0. 

It is easy to understand the reason why our modified algorithm did work : calling 

p — > R "position" 

(567) 

—x — > P "momentum" , 

the system can also be interpreted as an "ordinary" (i.e. massive and commutative) non- 
relativistic charged particle in the field of a Dirac monopole of strength 0, augmented with 
an inverse-square plus a Newtonian potential. This is the well-known "Mclntosh-Cisneros 
- Zwanziger" (MICZ) system [Mcintosh 19701 Zwanziger 196 8 1, for which the fine-tuned 



inverse-square potential is known to cancel the effect of the monopole, allowing for a 
Kepler- type dynamical symmetry [Mcintosh 1970} |Zwanziger 1968] . The angular momen- 



tum, (558), and the Runge-Lenz vector, (564), are, in particular, that of the MICZ problem 



[Mcintosh 19701 Zwanziger 1968| in "dual" momentum-space. 



The conserved quantities provide us with valuable information on the motion. Mimicking 
what is done for the MICZ case, we note that 

J-p=-0 (568) 

implies that the vector p moves on a cone of opening angle arccos (-9/J) . On the other 
hand, defining the conserved vector 

N = aJ-9K, (569) 

we construct the constant, 

N ■ p = 9(J 2 — 9 2 ) = 9L 2 , (570) 

so that the p-motion lies on the plane perpendicular to N. The trajectory in p-space belongs 
therefore to a conic section. 
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For the MICZ problem, this is the main result, but for us here our main interest lies in 



finding the real space trajectories, x(t). By (565), this amounts to find the [momentum-] 
"hodograph" of the MICZ problem. Curiously, while the hodograph of the Kepler problem 
is well-known, it is actually a circle or a circular arc, we could not find the corresponding 
result in the vast literature of MICZ system. 

Returning to our notations, we note that due to 



N-x = 0. 



(571) 



the x(t)-trajectories also belongs to an oblique plane, whose normal is N = a J — 6K . 
We can thus conveniently study the problem in an adapted coordinate system. One proves 
indeed that 



with A 2 = a 2 + 2UQ 2 , e 2 = a 2 + 2HJ 2 and L 2 = J 2 - 2 , 



(572) 



is a convenient orthonormal basis to study the x -trajectories. Here we recognize, in k, 
N/N in particular. 

• Firstly, projecting onto these axis, 



p z = p ■ k = 6L/\X\ = const 
^ p x =p-i, 

Py=P-3, 



(573) 



we find the equation 



Py + 



\e\at 

2\x\n 



i. 



A 2 /4-H 2 L 2 /2H ^ ^ 

which is the equation of a hyperbola or of an ellipse in momentum space, depending on the 
sign of H, positive or negative. For vanishing % one gets a parabola. This confirms what is 
known for the MICZ problem [Mcintosh 1970} Zwanziger 1968| , and is consistent with what 
we deduced geometrically. 



Next, projecting the x -motion onto the orthonormal basis (]572|) yields 
2\L\ 



X = x ■ i 



Y — x • j 
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A 


x ■ p 




6 


P 



Z — x ■ k — 



(575) 



An easy calculation leads to the equation 



a 2 L 2 



L 2 a 2 



-Y 



A 2 J 2 J 4 

which always describes an ellipse or an arc of ellipse, since 

x 2 = a 2 + me 2 > o . 



(576) 



(577) 



The center has been shifted along the axis % by the quantity ( — \e\Lj J 2 ) and the major axis 
is directed along j. Note that, unlike as in p -space, the x-trajectories are always bounded. 

When the energy is negative, % < 0, which is only possible when the Newtonian potential 
is attractive, a < , the x-trajectories are full ellipses. The origin is inside the ellipse : 



L= 1,H = -1/6, a = -1,8= 1 




-1.0 -0.5 0.0 0.5 1.0 



FIG. 8: % < and the Newtonian potential is attractive a < , so the trajectories describe a 
whole ellipse. 



When % > 0, which is the only possibility in the repulsive case a > , the origin is 
outside the ellipse so that only the right arc [denoted with the heavy line in the left side 
figure of ([9])] between the tangents drawn from the origin is obtained. However, positive 
hamiltonian H > , is also allowed for attractive Newtonian potential a < but in that 
event the origin is again outside the ellipse so that the x -trajectories are confined on the 
left arc of the ellipse [denoted with the heavy line in the right side figure of (|9j)] : 
For H = 0, the origin lies on the ellipse, and "motion" reduces to this single point : 
When the non-commutativity is turned off, 9 — > 0, the known circular hodographs of 
the dual Kepler problem are recovered. As a — > 0, the trajectory becomes unbounded, and 
follows the y-axis. 
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FIG. 9: The left side figure represents a right arc of an ellipse spanned by the x -trajectories for 
H > and a > 0. While the right side figure represents a left arc of an ellipse spanned by the 
x -trajectories for H > and a < 0. 



L= l,H = 0,g= -1,8= 1 




-1.5 -1.0 -0.5 0.0 0.5 



FIG. 10: x -trajectories degenerate to one single point for % = 0. 

So far, we only discussed classical mechanics. Quantization is now straightforward using 
the known group theoretical properties of the MICZ problem in dual space. The non- 
commutativity, alias monopole charge, 9 has to be an integer or half integer. This is indeed 
the first indication about the quantization of the non-commutative parameter. The wave 
functions should be chosen in the momentum representation, if)(jp). The angular momentum, 
J, and the rescaled Runge-Lenz vector, Kj y / 2~|W|, close into o(3, l)/o(4) depending on the 
sign of the energy. In the last case, the representation theory provides us with the discrete 
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energy spectrum, see (123), [in units h = 1} 

" 2 1 1 



E n 



a 
2^2' 



n = n r + - + (l+ 2 



1 + 



49 2 



(21 + 1) 



2 ' 



(578) 



where n = 0, 1, . . . , I = 0, 1, . . . , with degeneracy 



n 



e 2 = (n - e) (n + e). 



Note that the degeneracy always takes integer or half-integer value, as it should, since n 
and 9 are simultaneously integer or half-integer. The same result can plainly be derived 
directly by solving the Schrodinger equation in j>-space [Mcintosh 19701 Zwanziger 1968 



Also related to the MICZ system, calculation of energy levels of hydrogen atom using NC 
QED theory is discussed in [Chaichian 2001] . 

Moreover, the symmetry extends to the conformal o(4, 2) symmetry, due to the fact that 
the massless Poincare orbits with helicity 9 are in fact orbits of the conformal group, cf. 
[Cordani 1990] . 

Let us observe that in most approaches one studies the properties (like trajectories, sym- 
metries, etc.) of some given physical system. Here we followed the reverse direction: after 
positing the fundamental commutation relations, we were looking for potentials with re- 



markable properties. This leads us to the momentum-dependent potentials (563), realizing 



a Mclntosh-Cisneros-Zwanziger system [Mcintosh 1970j. Zwanziger 1968| in dual space. Un- 
like as in a constant electric field |Horvathy 2006| , the motions lie in an (oblique) plane. 
The particle is confined to bounded trajectories, namely to (arcs of) ellipses. 

The best way to figure our motions is to think of them as analogs of the circular 
hodographs of the Kepler problem to which they indeed reduce when the non-commutativity 
is turned off. For % < 0, for example, the dual motions are bound, and the velocity turns 
around the whole ellipse; for H > instead, the motion along a finite arc, starting from 
one extreme point and tending to the other one at the end of the arc, corresponds to the 
variation of the velocity in the course of a hyperbolic motion of a comet, or in Rutherford 
scattering, but in dual space. 



Our system, with monopole-type non-commutativity (550), has some remarkable proper- 
ties : 

Momentum-dependent potentials are widely used in nuclear physics, namely in the study 
of heavy ion collisions, where they correspond to non- local interactions [Gale 19871 TDas 20031 
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IDas 2004] . Remarkably, in non-commutative field theory, a 1/p 2 contribution to the propa- 
gator emerges from UV-IR mixing. 

The absence of a mass term should not be thought of as the system being massless; it is 
rather reminiscent of "Chern-Simons dynamics" [Dunne 1 990]. 

One can be puzzled how the system would look like in configuration space. Trying to 



eliminate the momentum from the phase-space equations (565) in the usual way, which 
amounts to deriving x with respect to time and using the equations for p, fails, however, 
owing to the presence of underived p in the resulting equation. This reflects the non-local 
character of the system. 

One can, instead, eliminate x using the same procedure, but in dual space. This yields 
in fact the equations of the motion of MICZ in dual momentum space, 

J 2 a ^ a - 7 /K««\ 

P=~7 + ^, P= -SP- ■ ( 579 ) 

p6 pZ pd pi 

Are these equations related to a theory with higher-order derivatives of the type 
[Lukierski 1997] ILukierski 2003j ? The answer is yes and no. The clue is that time is 
not a "good" parameter for Kepler- type problems, owing to the impossibility of expressing 
it from the Kepler equation |Cordani 2003] . This is also the reason for which we describe the 
shape of the trajectories, but we do not integrate the equations of the motion. A "better" 
parameter can be found along the lines indicated by Souriau |Souriau 1982] IBates 1989] and 



then, deriving with respect to the new parameter, transforms (579 ) into a fourth-order linear 
matrix differential equation, which can be solved. 

It is, however, not clear at all if these equations derive from some higher-order Lagrangian, 
and if they happen to do, what would be the physical meaning of the latter. 

The fourth-order equations do certainly not come from one of the type stated in Ref. 
[Lukierski 1997] ILukierski 2003j : the latter lives in fact in two space dimensions and has con- 
stant scalar non-commutativity 6, while our system is 3-dimensional and has a momentum- 



dependent vector 0(p), given in (550). 

It is tempting to ask if the relation to the "closest physical theory" with a momentum- 
dependent potential, namely nuclear physics, can be further developed and if similar (su- 
per) symmetries can be found also in nuclear physics. Once again, the answer seems to be neg- 
ative, though : while dynamical symmetries do play a role in nuclear physics |IachelloT9 93] , 
those used so far do not seem to be of a momentum-dependent Keplerian type. 
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VII. CONCLUDING REMARKS 



In this thesis, we developed a systematic method to search for hidden and (su- 
per) symmetries of several physics system. In some cases, like in the SUSY of the monopole, 
our recipe needs to be extended to fermionic degrees of freedom. In the case of the 
momentum-space monopole, we needed to adapt our recipe to the non-commutative struc- 
ture by interchanging the role of positions and momenta. In these models, as expected, the 
hidden symmetry of the Kepler-type is always related to the addition into the system of 
a fine-tuned inverse-square potential. This requirement appears clearly, in the van Holten 
algorithm, to be a consistency condition on the existence of a conserved Runge-Lenz-type 
vector. 

Having introduced the Abelian Dirac magnetic monopole-field; we studied, in particular, 
the classical geodesic motion of a particle in Kaluza-Klein-type monopole spaces and its 
generalization: the Gibbons-Hawking space. We derived the conditions under which the 
Killing tensors imply the existence of conserved quantities on the dimensionally reduced 
curved manifold. We observed that the Killing tensor generating the Runge-Lenz-type vec- 



tor, preserved by the geodesic motion, can be lifted to an extended manifold, namely, (210) 



and (216) [Duval 19 9 lj. As an illustration, we have treated, in detail, the generalized Taub- 
NUT metric, for which we derived the most general additional scalar potential so that the 
combined system admits a Runge-Lenz vector Gibbons 09/2006"] . Another example con- 
sidered is the multi-center metric where we have found a conserved Runge-Lenz-type scalar 



(207), in the special case of motions confined onto a particular 2-sphere. Moreover, from 
the Theorem III.4| we deduced, for N > 2, that no Runge-Lenz vector does exist in the 
case of A-center metrics. It is worth mentioning that apart from the generic importance of 
constructing constants of motion, namely in the confinement of particle to conic sections; 
the existence, in particular, of quadratic conserved quantity like Runge-Lenz vector yields 
the separability of the Hamilton- Jacobi equation for the generalized Taub-NUT metric and 
for the two-center metric. 

In the case of isospin-carrying particle in a non-Abelian Wu-Yang monopole field, we 
found the most general scalar potential such that the combined system admits a conserved 
Runge-Lenz vector. Indeed, it generalizes the fine-tuned inverse-square plus Coulomb poten- 
tial [Mcintosh 1970} Zwanziger 1968| , for a charged particle in the field of a Dirac monopole. 
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Following Feher, the result is interpreted as describing motion in the asymptotic field of a 
self-dual Prasad-Sommerfield monopole |Feher 19841 IFeher 1985| iFeher 1 986J. 

We also treated the case of the effective "truly" non-Abelian monopole-like field generates 
by nuclear motion in a diatomic molecule. This system is due to Moody, Shapere and Wilczek 



where despite the non- conservation of the electric charge (386), we surprisingly constructed, 
in addition to the "unusual" angular momentum (??), a new conserved charge (405). 

We remarked that Runge-Lenz-type vector plays a role also in SUSY. Indeed, we investi- 
gated the bosonic symmetries as well as the supersymmetries of a spinning particle coupled to 
a magnetic monopole field. The gyromagnetic ratio determines the type of (super) symmetry 



the system can admit : for the Pauli-like hamiltonian (412) M = 1 SUSY only arises for 
gyromagnetic ratio g = 2 and with no external potential, V = 0, confirming Spector's ob- 
servation |Spector| . We also derived additional supercharges, which are not square roots of 
the Hamiltonian of the system, though. A Runge-Lenz-type dynamical symmetry requires 
instead an anomalous gyromagnetic ratio, 



g = or g — 4 



with the additional bonus of an extra "spin" symmetry. These particular values of gyro-ratio 
come from the effective coupling of the form Fij =F e^D^, which add or cancel for self-dual 
fields, = eijkDk® |Feher 1988] . We found that the super- and the bosonic symmetry can 
be combined, but the price to pay is, however, to enlarge the fermionic space. This provides 
us with an A/" = 2 SUSY. 

We also applied the van Holten algorithm to a planar fermion in any planar magnetic 
field, i.e. one perpendicular to the plane. We shown, for ordinary gyromagnetic, that in 



addition to the usual supercharge (532) generating the supersymmetry, the system also 
admits another square root of the Pauli Hamiltonian |H orvathy 2005| happening due to the 
existence of a dual Killing tensor. 

A three-dimensional non- commutative oscillator with no mass term but with a cer- 
tain momentum-dependent potential is obtained when studying the hidden symmetry of 



a monopole-type non-commutativity |N gome 06/2010] . This oscillator system exhibits a 



conserved Runge-Lenz-type vector derived from the dual description in momentum space. 
The latter corresponds, but in dual space, to a Dirac monopole with a fine-tuned inverse- 
square plus Newtonian potential, introduced by Mcintosh, Cisneros, and by Zwanziger some 



126 



time ago. The resulting additional Kepler-type symmetry leads to the confinement of the 
particle's trajectories to bounded trajectories, namely to (arcs of) ellipses. When the non- 
commutativity is turned off, i.e. in the commutative limit, the motions reduce to the circular 
hodographs of the Kepler problem. It is worth mentioning that the momentum-dependent 
potentials which are rather unusual in high-energy physics, however, are widely used in 
nuclear physics, namely in the study of heavy ion collisions; they correspond to non-local 
interactions |Gale 1987] IDas 2003]. IDas 2004] . Moreover, in non-commutative field theory, it 
is remarkable that a 1/p 2 contribution to the propagator emerges from the UV-IR mixing. 
See in |Gubser 2001] . The absence of a mass term in the Hamiltonian describing this non- 
commutative oscillator should not be thought of as the system being massless; it is rather 
reminiscent of "Chern-Simons dynamics" [Dunne 1990j . 
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